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ABSTRACT 


Temperature  fluctuations  in  the  atmosphere  severely  limit  the  angular  resolution  of  earth 
bound  observation  facilities  to  around  1  arcsecond.  This  corresponds  to  an  effective,  coherent, 
af)erture  size  of  10  cm  even  though  the  telescope  may  have  a  2-4  m  primary  mirror. 
Understanding  the  spatial  and  temporal  distribution  of  atmospheric  optical  turbulence  is 
essential  to  maximize  the  performance  of  large  astronomical  telescopes.  This  thesis  made  use 
of  a  5  kHz  high  frequency,  short  range  Doppler  acoustic  sounder  to  investigate  the  first  100 
meters  of  the  mountain  boundary  layer  turbulence  above  the  Air  Force  Maui  Observation  Site, 
AMOS,  Haleakala,  Hi.  These  measurements  were  part  of  a  coordinated  site  evaluation  for  a 
proposed  4  m  telescope  to  be  built  at  AMOS  in  the  near  future. 

Tentative  results  revealed  significant  layering,  15-20  m  and  occasionally  thicker,  in  the 
turbulent  surface  layers  above  AMOS.  Add'tionally,  a  comparison  of  two  proposed  construction 
sites  near  the  top  of  Haleakala  showed  that  the  turbulent  surface  layer  tends  to  follow  the 
contours  of  the  mountain. 
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I.  INTRODUCTION 


A.  OPTICAL  TURBULENCE 

The  resolution  of  an  eistronomical  telescope  is  known  to  be  highly 
dependent  upon  localized  optical  turbulence  characteristics.  The  ability  to 
measure  eind  to  understand  such  phenomenon  can  enhance  the  optical  viewing 
quality  of  earth-bound  instruments.  One  method  of  monitoring  these  naturally 
occurring  effects  at  low  altitudes  (<200m),  is  through  the  use  of  am  acoustic 
echo-sounder.  A  variety  of  such  devices  have  been  proven  to  be  both  feasible 
and  accurate  [Refs.  1  through  4].  The  echo-sounder  used  to  obtain  the  data  in 
this  thesis  has  been  developed  over  the  past  4  years  by  Walters,  Moxcey, 
Weingartner,  Wroblewski,  and  McCrary  [Refs.  5  through  9]. 

Temperature  fluctuations  in  the  atmosphere  degrade  optical  seeing 
conditions  by  affecting  the  local  index  of  refraction.  The  scattering  of  acoustic 
energy  produced  by  an  echo-sounder  can  reveal  the  spatial  and  temporal 
distribution  of  these  fluctuations.  The  acoustic  echo-sounder  exploits  this 
scattering  phenomenon  to  model  optical  turbulence. 

B.  MISSION  OBJECTIVES 

The  mission  of  this  project  was  to  assist  in  the  site  assessment  for  a 
proposed  Air  Force,  4  meter  telescope  at  the  Air  Force  Maui  Obsen/ation  Station, 
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AMCo,  Haleakala,  Hi.  This  was  accomplished  in  two  phases.  The  first  of  these 
investigated  the  existence,  location,  and  significance  of  low  altitude  turbulence 
layers.  The  second  performed  a  quick-look  comparison  between  two  proposed 
construction  sites. 

Sounder  measurements  extended  between  26  June  and  2  July,  1991,  at 
both  proposed  sites.  Additionally,  the  1.6  meter  Itek  Compensated  Imaging 
System,  CIS,  and  the  1.2m  Hughes  seeing  monitor  gathered  simultaneous  stellar 
full  width  at  half  maximum  (FWHM)  data.  The  primary  site.  Site  A,  was  at  the 
summit  of  Haleakala  in  the  center  of  a  30  m  circular  ring,  15  m  higher  and  100 
m  to  the  south-east  of  the  current  facility.  The  secondary  site.  Site  G,  was  located 
35  m  to  the  south-east  of,  and  at  approximately  the  same  elevation  as,  the 
current  facility.  A  layout  of  the  AMOS  facility  is  given  in  Figure  1 . 
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Figure  1  AMOS  Facility  Layout  Showing  Locations  of  Proposed  4m 
Telescope  Sites. 
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II  BACKGROUND 


A.  ATMOSPHERIC  TURBULENCE 

An  optical  plane  wave  experiences  random  phase  perturbations  while 
passing  through  the  atmosphere.  The  cause  of  these  perturbations  stems  from 
the  presence  of  small  temperature  fluctuations  carried  by  turbulent  velocity  fields 
in  the  atmosphere.  These  phenomena  affect  the  local  index  of  refraction  directly, 
and  hence,  degrade  the  resolution  of  earth-bound  observation  facilities.  The 
following  two  parameters  arise  when  describing  this  degradation;  1)  the  refractive 
turbulence  structure  parameter,  and  2)  the  spatial  coherence  length  of  the 
atmosphere,  r^ . 

The  more  fundamental  of  these,  is  the  mean-square  statistical  average 
of  the  difference  between  the  indices  of  refraction  of  two  points.  Tatarski  [Ref.  10] 
includes  a  detailed  discussion  of  this  parameter  which  is  summarized  more 
recently  by  Walters  and  Kunkel  [Ref.  11]  as  follows, 

zli^  (1) 

where  n,  and  nj  are  the  indices  of  refraction  at  points  1  and  2  respectively,  and 
ri2  is  the  separation  between  points  1  and  2  in  the  atmosphere.  It  is  a  measure 
of  the  local  optical  turbulence  present  in  the  atmosphere  and  is  the  starting  point 
in  analyzing  that  turbulence.  It  is  important  to  note  that  the  presence  of  significant 
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optical  turbulence  does  not  imply  an  appreciable  degree  of  velocity  turbulence 
and  vice  versa  [Ref.9].  Regretfully,  is  a  difficult  parameter  to  measure  directly. 
It  is  easier  to  make  direct  measurements  of  another  atmospheric  structure 
parameter,  Cj^. 

The  thermal  structure  parameter,  Ct^,  is  the  mean  squared  temperature 
difference  between  two  points  in  space, 

(2) 

where  T,  and  Tj  are  the  temperatures  at  points  1  and  2  respectively.  As  will  be 
discussed  in  the  following  section,  an  acoustic  echo-sounder  measures  this 
parameter  indirectly.  Tatarski  [Ref.  10]  develops  a  relationship  between  C/  and 

'-'N  • 

Cw  =  (79xlO'®P/r^)2Cr  (3) 

where  P  is  the  atmospheric  pressure  in  millibars,  and  T  is  the  temperature  in 
kelvins.  Equation  (3)  ignores  a  contributing  factor  that  depends  upon  water  vapor 
concentration.  This  is  a  good  approximation  for  optical  wavelengths  under  dry, 
non-maritime  conditions  [Ref.  1 1  ]  such  as  those  present  at  AMOS. 

As  mentioned  previously,  the  spatial  coherence  length,  r„,  is  the  second 
optical  turbulence  parameter.  It  represents  the  magnitude  of  the  integrated. 


optical-turbulence  present.  It  depends  on  the  integral  of  over  the  distance  in 
question  [Ref.  11].  For  a  plane  wave, 

To  =  2.1  [  1.46ic2  C^(z)dz  ] -3/5  (4) 


where  k  is  the  wave  number  (2Tr/X),  and  L  is  the  optical  path  length  along  the 
vertical  directon.  Vertical  path  coherence  lengths  vary  widely  between  observation 
sites.  Historically  r^  values  can  range  from  a  few  centimeters  for  a  poor  site  to 
approximately  30  cm  for  a  world  class  site.  The  following  sections  present  data 
in  the  form  of  coherence  lengths  integrated  over  several  different  optical  path 
lengths.  In  comparing  r^’s  between  sites  it  is  important  to  realize  that  the 
atmospheric  turbulence  above  the  site  is  not  the  only  variable  that  degrades 
seeing  ability.  Of  equal  importance  are  the  artificially  induced  conditions 
surrounding  and  within  a  facility  [Ref.  12].  Some  of  the  more  prevalent  of  these 
sources  are  ventilation  systems,  electrical  equipment,  and  convective  heat  from 
surrounding  buildings,  tarmacs,  and  optical  components. 

B.  ECHO-SOUNDER  OPERATION 

Acoustic  echo-sounders  use  a  transmitted  pulse  of  acoustic  energy  to 
investigate  the  properties  of  the  atmosphere  in  much  the  same  way  as  a  ship- 
borne  sonar  detects  the  presence  of  a  submarine.  The  temperature  structure  of 
the  atmosphere  backscatters  a  portion  of  the  transmitted  energy  that  the  echo- 
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sounder  receives.  The  echosonde  equation,  also  known  as  the  radar  equation, 
relates  the  scattered  acoustic  signal  to  a  number  of  atmospheric  and  instrumental 
parameters.  This  relationship,  based  on  work  by  Tatarski  [Ref.  10],  and  discussed 
by  Neff  [Ref.  2]  for  adaptation  to  an  echo-sounder  is, 

Pg  =  P-rEj^Ej.  (ct/2)  {a  ^{R,  f))  {aG/R^)  (5) 


where, 

•  Pr  is  the  acoustic  power  returned  form  a  range  R. 

•  I  ,  is  the  acoustic  power  transmitted  at  frequency  f. 

•  Er  is  the  efficiency  in  converting  from  acoustic  to  electrical  power. 

•  Ej  is  the  efficiency  in  converting  from  electrical  to  acoustic  power. 

•  T  is  the  length  of  the  acoustic  pulse  (sec). 

•  R  is  the  range  to  the  scattering  volume  (meters). 

•  c  is  the  local  speed  of  sound  (m/s). 

•  A  is  the  area  of  the  antenna  (m^). 

•  G  is  the  effective  aperture  factor  for  the  antenna. 

•  exp(-2aR)  is  the  round  trip  power  loss,  where  a  is  the  average  attenuation 
factor  over  the  path. 

•  a(f)  is  the  scattering  cross  section  per  unit  volume. 

As  mentioned  in  the  previous  section,  the  acoustic  sounder  measures 
atmospheric  turbulence,  Cj^,  values  from  the  scattering  cross  section,  a.  Tatarski 
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[Ref.  1 0]  shows  the  backscattering  cross  section  is  proportional  to  Cy®  through 
the  following  relationship, 


a(R,f)  =  0.0039  (Cr/T^)  (®) 

where, 

•  T  is  the  absolute  temperature  in  kelvins, 

•  k  is  the  acoustic  wave  number  =  2n/k, 

•  X  is  the  acoustic  wavelength. 

This  is  the  case  only  for  backscattered  acoustic  energy,  that  is,  energy  scattered 
through  an  angle  of  tt.  For  other  angles,  the  amount  of  energy  scattered  depends 
upon  the  velocity  structure  parameter  of  the  atmosphere  [Ref.  10].  Combining 
Eqs.  (5)  and  (6),  and  solving  for  C/  gives, 


^2  _  / _ _ \  /  Pr  \p2  g(2oR) 

^  0.0039Pj£j.^^G  /c^/3  (Ct)  (Pj.) 


(7) 


Several  of  the  factors  appearing  Equation  (7),  such  as  Et,  Er,  G  and  a,  are 
difficult  to  determine  accurately.  Therefore,  a  certain  degree  of  error  exists  with 
an  acoustic  sounder  that  must  be  dealt  with  when  calibrating  such  a  device.  The 
following  section  discusses  these  effects. 

1.  Acoustic  Sounder  Calibration 

The  advantages  of  an  echo-sounder  lie  in  its  simplicity  of  operation, 
portability,  and  its  ability  to  map  the  spatial  distribution  of  turbulence.  Regretfully, 
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calibration  of  the  data  from  an  acoustic  sounder  is  complex.  A  more  common 
means  of  measuring  temperature  variations  as  a  function  of  altitude  is  through 
the  use  of  high  speed  temperature  probes.  Mounting  these  probes  at  specific 
heights  along  a  rigid  tower  can  provide  direct  temperature  measurements. 
However,  such  devices  lack  the  portability  of  the  echo-sounder  and  are  costly  to 
build  and  to  maintain. 

Historically,  acoustic  echo-sounders  are  considered  to  be  semi- 
quantitative  devices.  Inherent  uncertainties  include  such  factors  as  difficulties  in 
measuring  transducer  efficiencies,  determination  of  the  effective  aperture  factor, 
G,  and  in  obtaining  a  precise  attenuation  factor,  a,  over  the  optical  path  length 
in  question.  As  with  any  experimental  system,  judicious  selection  of  equipment, 
close  monitoring  of  operating  parameters,  and  strict  experimental  procedures  will 
minimize  these  effects. 

Figure  2  is  a  plot  illustrating  the  dependence  of  the  attenuation  factor 
on  both  water  vapor  pressure  and  frequency.  Determining  a  precise  value  for  the 
atmospheric  attenuation  factor  as  a  function  of  altitude  is  extremely  complicated. 
Tower  mounted  humidity  probes  could  measure  the  height  dependence  but  then 
the  portability  benefit  of  the  sounder  is  lost.  In  order  to  solve  this  problem  an 
approximation  must  be  made.  One  such  method  is  to  take  a  ground  level  reading 
and  to  assume  it  to  be  the  average  value  over  the  path  length.  Since  this  project 
was  concerned  with  low-altitude  measurements,  this  approach  was  taken. 
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Figure  2.  Attenuation  vs  Water-Vapor  Pressure  for  Temperature  = 
20°C  and  Atmospheric  Pressure  =  1013  mb. 


The  inner  and  outer  scales  determine  the  dimensions  of  the  thermal 
turbulence.  These  refer  to  the  smallest  and  largest  scale  sizes  of  the  turbulence 
structure,  respectively.  Typical  values  for  these  parameters  range  from  millimeters 
to  hundreds  of  meters.  These  happen  to  be  of  the  same  order  of  magnitude  as 
acoustic  wavelengths.  Acoustic  sounders  measure  optical  turbulence  effectively 
because  of  this  match  between  acoustic  wavelength  and  characteristic  scale  size, 
combined  with  the  fact  that  the  acoustic  scattering  cross  section  is  orders  of 
magnitude  larger  than  the  optical  scattering  cross  section. 
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As  mentioned,  one  of  the  more  common  means  of  measuring  thermal 
turbulence  profiles  is  through  the  use  of  tower  mounted  temperature  probes.  An 
acoustic  sounder  may  be  calibrated  by  positioning  it  in  the  vicinity  of  such  a 
tower  and  comparing  simultaneous  measurements.  In  doing  so,  the  difference 
in  scale  sizes  measured  by  each  device  will  introduce  some  degree  of 
discrepancy  between  the  two  sets  of  readings. 

An  acoustic  echo-sounder,  unlike  a  temperature  probe  system,  is 
susceptible  to  inner  scale  errors  due  to  the  finite  dimensions  and  wavelength 
spectrum  of  its  acoustic  pulse.  Tatarski  [Ref.  10]  shows  that  the  Bragg  scattering 
expression  nk  =  2d  sin0  determines  the  turbulent  scale  size  probed  by  an 
acoustic  sounder.  For  the  5  kHz  acoustic  sounder  used  here,  this  scale 
corresponds  to  3.4  cm,  which  is  between  the  inner  and  outer  scales  of 
turbulence.  As  an  illustrative  example  consider  the  set  of  conditions  such  that  the 
mean  wind  speed  is  6  m/sec.  The  corresponding  inner  scale  length,  1^,,  of  the 
turbulence  present  for  these  winds  was  calculated  by  Ochs  and  Hill  to  be  «  3 
mm.  [Ref.  13]  using, 

ic  =[-]"/*  (8) 

€ 

where  v  is  the  kinematic  viscosity  and  €  is  the  turbulent  dissipation  rate.  Figure 
3  is  a  plot  of  the  spatial  power  spectrum  of  temperature  fluctuations,  (Dy,  versus 
scaled  wave  number,  kt]  taken  from  Reference  14  (both  are  dimensionless).  The 
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dotted  line  represents  the  Kolmogorov  spectrum  including  dissipation  for  velocity 
fluctuations.  The  solid  line  is  the  theoretical  model  of  Hill  [Ref.  14]  and 
corresponds  to  the  spectrum  of  temperature  fluctuations.  The  most  significant 
feature  of  this  plot  is  the  bump  present  in  the  Hill  model.  This  occurs  at  the  edge 
of  the  inner  scale  of  turbulence  and  arises  from  residual  temperature  gradients 
as  the  turbulence  enters  the  viscous  convective  region.  Viscosity  dissipates  the 
kinetic  energy  of  the  turbulence  while  temperature  gradients  diffuse  at  a  slower 
rate.  The  scaled  wave  number  consists  of  k,  which  is  equal  to  27r/scale  length, 
and  rj,  which  is  equal  to  1^/7. 14  (for  air).  Using  the  sounder  wave  length  of  3.4 
cm,  the  calculated  scaled  wave  number  is  approximately  7.8  x  10^.  From  Figure 
3  it  is  evident  that  this  point  lies  along  the  viscous-convective  bump  and 
corresponds  to  a  O^.  of  »  .88.  This  translates  to  an  enhancement  of  ~  20% 
compared  to  the  other  model.  The  position  of  this  bump  depends  on  the  inner 
scale  which  is  sensitive  to  the  local  dissipation  rate  €.  From  this  error  source 
alone,  C/  values  measured  by  the  acoustic  sounder  at  5kHz  should  be 
approximately  20%  greater  than  those  recorded  by  a  temperature  probe.  Of 
course,  this  calculation  was  conducted  only  for  a  specific  inner  scale.  This  value 
of  3  mm  was  chosen  since  it  was  representative  of  the  conditions  encountered 
at  AMOS. 
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Figure  3.  Spatial  Power  Spectrum,  Oj,  of  Temperature  Fluctuations  in  Air 
(both  axes  are  dimensionless).  The  Solid  Curve  is  the  Accurate  Model;  the 
Dashed  Cun/e  is  the  Tatarski  Model,  (from  Ochs  and  Hill  [Ref.  13]) 
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III.  EQUIPMENT 


A.  OVERVIEW 

The  majority  of  the  components  of  the  echo-sounder  used  in  this  project 
were  identical  to  those  described  by  Moxcey  [Ref.  6]  and  McCrary  [Ref.  9].  The 
system  as  a  whole  has  been  used  by  Walters  over  the  past  4  years  for  a  variety 
of  projects.  The  hardware  has  undergone  minor  refinements  during  this  time, 
while  the  software  has  undergone  three  major  revisions.  Figure  4  is  a  schematic 
of  the  entire  echo-sounder  system. 

The  acoustic  heart  of  the  system  is  a  19  element  hexagonally  shaped  array 
of  piezo-electric  speakers.  Figure  5  is  a  sketch  of  this  array.  The  individual 
elements  are  Motorola  KSN  1005A  speakers  and  were  chosen  for  their  resonance 
properties  as  well  as  their  high  efficiencies  in  both  the  transmit  and  receive 
modes.  The  array  beam  lobe  pattern  was  studied  extensively  in  an  anechoic 
chamber  and  found  to  be  highly  directional  [Ref.  6].  The  experimentally 
determined  parameters  associated  with  this  transducer  antenna  array  are: 

•  efficiencies:  E-r  =  Er  =  0.5, 

•  resonant  frequency  =  5kHz, 

•  area  of  antenna  =  A  =  0.0866  m^, 

•  effective  aperture  factor  =  G  =  .4. 


14 


To  accomplish  the  project’s  mission,  the  sounder  had  to  operate  in  a  mode 
that  would  provide  high-resolution  data  over  the  0  to  100  meter  altitude  range. 
From  Figure  2,  it  is  apparent  that  the  attenuation  of  an  acoustic  signal  by  the 
atmosphere  increases  with  frequency.  This  suggests  using  a  low  operating 
frequency.  However,  the  resolution  of  the  sounder  depends  on  wavelength  size 
and  is  inversely  proportional  to  the  frequency.  Additionally,  driving  the 
transducers  at  their  resonant  frequency  of  5  kHz  provides  a  high  efficiency  of 
nearly  50%.  However,  exploiting  this  resonance  condition  also  has  a  slight 
drawback.  After  transmitting  the  5  kHz  pulse,  the  transducers  continue  to  ring 
momentarily  as  does  any  oscillator  when  driven  at  resonance.  The  effect  of  this 
ringing  was  to  create  a  brief  post-transaction  period  when  reception  was 
impossible.  The  consequence  was  a  minimum  range  of  »  5-7  meters.  The 
following  acoustic  pulse  parameters  were  used  in  analyzing  the  7.5  to  100  meter 
layer  at  AMOS. 

•  frequency  of  sinusoidal  pulse  =  5kHz, 

•  Pulse  length  =  75  cycles, 

•  Magnitude  of  Pulse  =  40  volts, 

•  Resolution  =  (c*t)/2  =  (.5)(75*340m/s)/(5kHz)  =  2.55  meters. 

Past  experience  has  shown  that  this  is  an  optimum  choice  of  operating 
parameters. 
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B.  RECENT  UPGRADES 


Upgrades  have  been  made  in  the  areas  of  computing  power,  mass  data 
storage,  Doppler  processing  and  acoustic  attenuation.  The  control  of  all 
components  is  now  handled  by  a  Hewlett-Packard  9836C  P.  C.,  vice  the  previous 
H-P  217.  This  is  an  upgrade  over  the  previous  system  in  that  speed  of  operation 
has  been  significantly  increased  through  the  addition  of  a  Newport  Digital  33MHz 
68020  control  processing  board.  The  standard  H.P.  3.5"  disk  drive  has  been 
replaced  with  a  Bering  600  megabyte  optical  disk  drive.  This  enables  storage  of 
individual  pulses  of  raw  turbulence  data  which  were  previously  recorded  only  as 
line  printer  plots.  The  massive  storage  capsdDility  of  the  optical  disks  provides  the 
ability  to  record  several  weeks  of  data  on  the  same  disk. 

The  last  equipment  modification  is  in  the  preamplifier.  This  change  was 
made  in  order  to  avoid  R.F.  interference  problems  that  were  experienced  earlier 
[Ref.  6].  The  new  preamp  can  be  operated  at  several  gains  (measured  to  be  987, 
1906,  or  3895)  and  was  designed  by  D.  L.  Walters  and  built  by  D.  Glarowicz.  For 
the  data  presented  in  this  paper,  a  gain  of  1 906  was  used  throughout. 

The  computer  software  provides  instrument  control  and  real  time  data 
processing.  Acoustic  pulse  length,  number  of  cycles,  pulse  strength,  filtering  of 
obvious  noise  signals,  data  storage  and  data  display  are  all  functions  controlled 
by  the  sounder  data  collection  program.  This  software,  [Appendix  A],  has 
undergone  numerous  revisions  since  the  first  version  was  drafted  in  1987  by 
Wroblewski  and  Walters  [Ref.  8].  Of  these  changes,  three  stand  out  as  being 
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especially  significant.  The  first  of  these  involves  the  addition  of  a  Doppler 
algorithm.  This  was  written  by  D.  L.  Walters  and  investigated  by  McCrary  [Ref.  9]. 
It  provides  for  the  collection  of  horizontal  wind  speeds  by  tilting  the  sounder 
slightly  off  the  vertical  and  measuring  the  Doppler  shift  in  the  backscattered 
signal.  This  enables  additional  turbulence  information  regarding  the  velocity  field 
of  the  atmosphere  to  be  recorded.  The  second  significant  software  change 
involves  the  gain  of  the  reflected  signal.  Previously,  a  constant  gain  of 
approximately  100,000  was  used  to  amplify  the  reflected  signal  over  the  entire 
range  of  altitude  concerned.  This  often  caused  some  overloading  on  the  filter 
during  periods  of  excessive,  low-altitude  turbulence,  but  was  necessary  in  order 
to  obtain  measureable  return  from  heights  of  100  m.  Now  there  is  a  gain  change 
accomplished  during  acquisition.  From  0  to  30  meters,  the  Infotek  AD-200  analog 
to  digital  converter  provides  unity  gain  to  the  signal  received  from  the  filter.  From 
30  to  100m,  the  AD-200  imparts  an  additional  gain  multiple  of  10.  This  extends 
the  operating  range  without  causing  an  overload  at  short  ranges.  Numerically, 
this  translates  to  the  following  gains.  Over  the  altitude  range  of  0  to  30  meters, 
the  total  gain  =  (pre-amp  gain)  (filter  gain)  (soft  ware  gain)  =  (1906)(10)(1)  = 
19,060.  Over  the  altitude  range  from  30  to  100  meters,  the  total  signal  gain  is  = 
(1906)  (10)  (10)  =  190,600.  A  16  bit  A-D  converter  would  alleviate  the  need  for 
such  a  gain  change,  but  at  the  time  of  this  project  none  was  readily  available. 
The  third  of  these  software  refinements  involves  a  modification  to  the  water  vapor 
acoustic  attenuation  algorithm,  based  on  a  recent  publication  by  H.  Bass  [Ref. 
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15].  The  effect  of  this  modification  is  to  significantly  change  (by  as  much  as  a 
factor  of  2)  the  absorption  coefficient  at  very  low  and  very  high  humidity  values 
as  compared  to  the  previous  algorithm.  In  the  middle  humidity  range,  there  is 
little  change  from  the  previously  determined  value  [Ref.  2]. 
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A.  CALIBRATION 


1.  Equipment 

Remote  sensors  must  be  calibrated  to  provide  accurate  data.  In  the 
case  of  acoustic  echo-sounders,  this  is  usually  accomplished  through  the  use  of 
temperature  probes  positioned  along  a  tower  at  specified  heights.  This  provides 
a  means  for  obtaining  an  independent  C/  vs.  altitude  profile  for  comparison. 
Ideally,  probes  are  positioned  along  the  tower  over  the  entire  altitude  range  of  the 
sounder  being  calibrated.  However,  due  to  cost,  equipment,  and  time  constraints, 
the  calibration  conducted  for  this  project  consisted  of  a  single  point  comparison 
at  heights  of  10  and  14  meters.  These  heights  were  chosen  for  two  reasons.  First, 
they  were  consistent  with  reasonable  pedestal  heights  for  the  AMOS  4m  facility. 
Secondly,  the  lower  altitude  region  was  of  special  interest  since  it  allowed 
comparisons  with  optical  data  from  the  1.6m  CIS  and  1.2m  Hughes  telescopes. 
This  calibration  experiment  was  conducted  on  the  roof  of  Spanagel  Hall,  Naval 
Postgraduate  School,  Monterey^  California.  Figure  6  shows  the  equipment 
positioning  for  the  calibration.  To  conduct  an  accurate  calibration,  the  sounder 
and  the  temperature  probe  should  be  measuring  the  identical  atmospheric 
volume.  However,  since  the  probe  itself  will  scatter  acoustic  energy,  it  must  be 
offset  from  the  sounder.  As  depicted  in  Figure  (6),  the  predominant  wind  direction 


Figure  6.  Equipment  Positioning  for  Acoustic  Sounder  Calibration. 
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during  the  calibration  would  appear  to  allow  both  devices  to  measure  the  same 
Cj‘,  while  minimizing  building  or  probe  induced  turbulence  contributions. 

The  differential  thermocouple  probe  used  is  one  which  was  designed 
by  Walters  and  investigated  by  Olmstead  [Ref.  15].  It  consists  of  two  0.00254  cm. 
diameter  Copper-Constantan  thermocouples  rigidly  attached  to  the  ends  of  a  1 
meter  aluminum  rod.  This  comprises  a  device  that  is  capable  of  determining  a 
Ct^  value  within  a  sensitivity  of  =  10"*  (i^/m^®).  Figure  7  contains  a  baseline 
measurement  of  this  sensitivity  taken  in  the  Naval  Postgraduate  School  anechoic 
chamber. 

2.  Findings 

Data  were  taken  on  three  days,  Oct.  1 0,  23,  and  24,  when  wind  and 
weather  conditions  were  appropriate.  Initial  data  taken  at  10  m  indicated  a  major 
discrepancy  between  the  C/  values  measured  by  the  the  sounder  and  the  probe. 
As  discussed  in  the  Background  section,  it  was  anticipated  that  there  would  be 
a  15-20%  difference  between  the  two  values  due  to  scale  size  effects.  However, 
comparing  the  two  data  sets  revealed  that  the  sounder  data  was  a  factor  of  2-4 
times  smaller  than  the  thermal  probe  results.  Additionally,  there  were  several 
other  trends  in  the  data  which  were  not  understood.  Figures  8,  9,  and  1 0  contain 
these  results.  Figure  8  contains  the  only  set  of  data  taken  at  a  height  of  10 
meters.  Notice  that  the  data  points  in  this  plot  appear  to  be  in  closer  agreement 
than  those  of  the  other  two.  The  purpose  for  taking  measurements  at  two 
altitudes  was  to  ascertain  some  measure  of  the  dynamic  range  of  the  sounder. 


Ct^Z  (lag  average) 


B  Oct  1391  BflSELIflEOl 


T  1  me  »'  UT  C  > 


Figure  7.  Differential  Thermocouple  Probe  Sensitivity  Baseline 
Sensitivity  Measurerr  ent. 


Local  Time 

Figure  8.  Acoustic  Sounder  Calibration  Data  for  10  October,  1991  at 
a  Height  of  10m. 
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Local  Time 

Figure  9.  Acoustic  Sounder  Calibration  Data  for  23  October,  1991  at 
a  Height  of  15m. 


Figure  10.  Acoustic  Sounder  Calibration  Data  for  24  October,  1991  at 
a  Height  of  15m. 
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The  fact  that  a  noticeable  difference  existed  over  such  a  small  range  of  altitudes 
was  intriguing.  Initial  investigation  into  the  cause  of  these  discrepancies  centered 
around  possible  hardware  causes.  This  included  replacement  of  both 
thermocouples,  detailed  measurements  of  both  the  sounder  and  probe  amplifier 
parameters,  and  replacement  of  several  lower  efficiency  transducer  elements. 
Discovering  no  significant  changes,  the  investigation  shifted  towards  possible 
contamination  sources.  Was  there  some  contamination  of  the  measurements 
from  building  induced  effects  at  the  lower  altitude?  Upon  reviewing  the 
equipment  positioning,  and  obsen/ing  the  trends  in  the  prevailing  winds,  the 
following  conclusions  were  made. 

To  record  data  at  a  height  of  10m  the  thermocouple  probe  was 
positioned  over  the  sounder  at  nearly  the  same  height  as  the  upper  roof  line  in 
Figure  6.  As  discussed,  the  set  up  was  chosen  so  that  the  prevailing  winds  would 
allow  both  devices  to  measure  the  same  uncontaminated  atmospheric  volume. 
However,  upon  closer  analysis  the  10  m  location  was  found  to  be  sensitive  to 
eddy  formation  near  the  roof.  This  artificially  induced  turbulence  would 
undoubtedly  cause  the  sounder  to  measure  different  Cj®  values  at  10m  than  at 
15m  when  compared  to  the  probe.  This  reasoning  explained  the  discrepancy 
between  the  10m  and  15m  sounder  data  but  still  didn’t  explain  the  persistent 
factor  »  4  between  the  two  devices  In  the  15m  data.  Additionly,  the  23  and  24 
Oct.  data  still  contained  other  inconsistencies  that  required  explanation. 
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Both  days  contained  differences  between  successive  sets  of  data  that 
varied  by  a  factor  of  two  or  more.  To  explain  these,  some  interaction  between  the 
building  and  the  wind-borne  turbulence  structure  was  suspected.  The  observed 
weather  conditions  for  each  day  were  similar;  clear  and  sunny  in  the  morning, 
gusty  afternoon  winds  becoming  more  directional  as  the  afternoon  passed,  and 
becoming  near  calm  near  sunset  (approx.  17:30  local  time).  Accurate  calibration 
required  that  the  winds  come  consistently  from  the  ocean.  Searching  for  trends 
in  the  data  revealed  that  the  measurements  differed  by  a  consistent  factor  of  four 
during  several  hours  of  the  day.  Upon  reviewing  the  weather  observations,  the 
winds  blew  steadily  from  the  ocean  during  these  hours.  Prior  to  these  periods  the 
winds  were  gusty.  Following  these  periods  the  winds  became  light  and  variable. 
Both  of  these  conditions  would  influence  the  turbulence  detected  by  both 
devices.  The  atmospheric  volume  each  was  measuring  was  not  only  different  but 
was  also  artificially  contaminated.  Eddying  effects  as  well  as  convective 
interaction  with  the  rooftop  and  exterior  building  walls  during  these  periods 
invalidated  the  measurements.  Based  on  these  arguements,  the  only  data  from 
the  two  days  considered  trustworthy  were  those  taken  during  the  mid-afternoon 
hours  when  sea  breeze  wind  conditions  prevailed.  Subsequently,  a  detailed 
review  of  the  sounder  operating  software  was  conducted.  The  intent  of  this  was 
to  find  an  explanation  for  the  consistent  factor  of  4  ratio  between  the  two  devices 
during  these  valid  periods. 


In  programmimg  Equation  7  into  the  data  processing  code,  an  error 
was  discovered  in  the  placement  of  the  receive  and  transmit  efficiency  factors. 
This  mistake  lowered  the  measured  sounder  values  by  a  factor  of  four.  By 
comparing  the  sounder  and  probe  calibration  data,  it  was  found  that  the  sounder 
measured  a  Cj®  value  «  3.8  times  lower  than  the  probe.  This  factor  of  3.8  is 
consistent  with  the  factor  of  four  coding  error  after  including  the  correction  for  the 
inner  scale  bump  in  Figure  3.  The  3.8  value  was  used  as  a  calibration  for  all  the 
AMOS  data,  and  is  appropriate  for  all  data  previously  recorded  by  this  instrument. 
Figures  11  and  12  show  the  corrected  plots  of  the  calibration  data  from  23  and 
24  Oct.  1991. 


Figure  11.  Acoustic  Sounder  Calibration  Data  for  23  Oct.,  1991 
Corrected  Using  Factor  of  3.8  x  Sounder  Data. 
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Figure  12.  Acoustic  Sounder  Calibration  Data  for  24  Oct.,  1991 
Corrected  Using  a  Factor  of  3.8  x  Sounder  Data. 

B.  AMOS  FIELD  MEASUREMENTS 

1.  Turbulence  Profile 

The  AMOS  mission  required  a  high-resolution,  low-altitude  turbulence 
profile.  Hence,  the  sounder  operating  parameters  were  chosen  to  collect  accurate 
turbulence  readings  over  a  range  of  7.5  to  100  meters.  As  discussed  in  the 
Equipment  chapter,  the  resolution  of  the  sounder  in  this  mode  was  **  2.6  meters, 
which  was  considered  ideal  for  the  ranges  in  question.  Some  representative 
examples  of  the  data  taken  at  Site  A  are  given  in  Figures  13  through  16.  The 
uppermost  of  the  plots  displayed  in  each  set  is  a  time  vs.  altitude  plot  of  the 
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backscattered  signal  received  by  the  sounder.  The  intensity  of  the  signal  received 
is  proportional  to  the  degree  of  shading  of  the  plot.  The  lower  of  the  two  plots  in 
each  example  is  a  15  min.  time-averaged  plot  of  the  measured  value  vs. 
altitude.  Most  of  the  data  was  for  Site  A,  the  primary  proposed  site  under 
consideration.  In  comparing  these  four  examples,  it  is  evident  that  the  conditions 
present  varied  considerably  from  day  to  day.  The  profiles  illustrated  in  Figures  13 
and  14  are  representative  of  the  typical  conditions  observed  at  AMOS.  Though 
Figures  15  and  16  are  atypical,  it  is  interesting  to  note  that  even  during  these 
high  turbulence  periods  there  is  little  if  any  signal  return  above  90  meters. 
Appendix  B  contains  additional  examples  of  sounder  data  plots. 

Moving  the  acoustic  sounder  antenna  between  sites  A  and  G  on  6 
different  occasions  provided  a  rough  intercomparison  between  the  two  locations. 
This  was  done  to  reveal  any  possible  gross  differences  between  the  two  sites. 
Figures  1 7  and  1 8  contain  a  few  representative  examples  of  these  Site  G  plots. 
No  significant  differences  appeared  between  the  two  locations  even  though  Site 
G  was  located  1 0  m  lower  than  Site  A  and  on  the  back  side  of  the  mountain 
peak.  Although  this  was  unexpected,  it  was  consistent  with  observations  of  fog 
flowing  over  the  mountain.  The  fog,  and  apparently  the  turbulent  boundary  layer, 
followed  the  contour  of  the  mountain  for  moderate,  6-10  m/s  wind  speeds. 

2.  vs.  Altitude 

As  well  as  measuring  the  turbulent  boundary  layer  profile,  an  equally 
important  question  was  finding  the  relative  contribution  of  the  first  1 00  m  to  the 
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Figure  14.  Sample  #2  of  Acoustic  Sounder  Data  Plots. 
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Figure  15.  Sample  #3  of  Acoustic  Sounder  Data  Plots. 
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Figure  16.  Sample  #4  of  Acoustic  Sounder  Data  Plots. 
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Figure  17.  Example  of  Site  G  Acoustic  Sounder  Data  Plot. 


entire  integral.  This  information  reveals  the  advantages  of  placing  the  4  m 
instrument  on  a  tower  or  pedestal  and  allows  a  tradeoff  of  optical  quality  vs.  cost 
of  a  tower.  Both  the  Kitt  Peak  and  Mauna  Kea  telescopes  are  4  m  instruments 
with  towers,  50  m  and  20  m  high  respectively,  that  demonstrate  the  feasibility  of 
this  option. 

Breaking  Equation  4  into  two  separate  integrals  over  the  optical  path 
length  gives  a  relation  for  determining  spatial  coherence  length,  r^,  as  a  function 
of  height,  z,  in  the  low  altitude  layer, 

-  H*  [  f^^°C^(z)dz+  f  C^j(z)  dz  ] (9) 

J7.5  JlOO 

where 

•  z  is  the  height  above  ground, 

•  is  the  measured  structure  parameter, 

•  is  the  structure  parameter  for  the  atmosphere  above  100  meters, 

•  H  is  the  product  of  all  constants  from  Equation  4  for  an  optical  wavelength 
of  500nm  and  is  =  5.0647  x  10®  (m'®'®). 

Equation  9  requires  an  approximation  regarding  the  turbulence 
contribution  from  the  second  integral.  The  Canada-France-Hawaii  Telescope, 
CFHT,  located  on  Mauna  Kea  provides  historical  data  to  achieve  this  [Ref.  12]. 
Because  CFHT  is  geographically  close  to  AMOS,  this  data  is  representative  of  the 
seeing  conditions  above  4  km.  Since  CFHT  averages  0.8  arcsec  FWHM  including 
its  surface  boundary  layer,  mirror  and  dome  effects,  a  20  cm  coherence  length 
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is  reasonable  for  AMOS  starting  at  100  m  above  the  surface.  Figures  19  through 
22  represent  plots  constructed  using  a  20  cm  coherence  length  in  Equation  (9) 
for  four  different  pedestal  heights  \within  the  7.5  to  100  m  layer.  Because  of  a 
combination  of  equipment  malfunctions  and  rainy  weather,  reliable  sounder  data 
is  missing  for  the  evenings  of  25,  27,  and  30  June. 

3.  Comparison  to  Optical  Data 
a.  AMOS  Support 

During  the  periods  of  sounder  operation,  both  the  1.6m  CIS  and 
the  1 .2  m  Hughes  seeing  monitor  at  AMOS  collected  optical  data  for  comparison. 
These  data  were  in  the  form  of  stellar  image  full  width  at  half  maxima,  FWHM, 
and  provides  a  means  for  confirming  the  sounder  results.  Comparing  the  optical 
and  acoustic  data  should  reveal  the  presence  of  artificially  induced  turbulence 
within  the  1.6  and  1.2  m  facilities.  Reference  12  is  a  recent  paper  discussing  the 
effects  of  such  turbulence  at  the  CFHT.  They  discovered  a  significant  amount  of 
degradation  from  such  factors  as  ventilation  systems,  wind  turbulence  across  the 
dome,  or  inhomogeneous  heating/cooling  of  the  mirror  itself.  The  following 
chapter  addresses  comments  regarding  the  effects  of  such  factors  on  the  1 .6m 
device.  To  compare  the  optical  and  sounder  results,  the  optical  and  acoustic  data 
must  be  consistent. 
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Figure  19.  Acoustic  Sounder  Data  for  the  Evening  of  26-27  June  1991 
Showing  the  Turbulence  Layering  Structure  at  AMOS. 
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Figure  20.  Acoustic  Sounder  Data  for  the  Evening  of  28-29  June,  1991 
Showing  the  Turbulence  Layering  Structure  at  AMOS. 
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Figure  21.  Acoustic  Sounder  Data  for  the  Evening  of  29-30  June,  1991 
Showing  the  Turbulence  Layering  Structure  at  AMOS. 
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Figure  22.  Acoustic  Sounder  Data  for  the  Evening  of  01-02  July  1991 
Shovnng  the  Turbulence  Layering  Structure  at  AMOS. 
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b.  Zenith  Angle  Correction 


The  first  step  in  correcting  the  optical  data  for  comparison  is  a 
zenith  angle  correclion.  This  is  a  simple  trigonometric  technique  that  involves 
division  of  the  optical  measurement  by  (seed)®'*  (where  6  is  the  respective  zenith 
angle)  to  rescale  the  data  for  the  slant  path  length  through  the  atmosphere. 
Caution  must  be  exercised  when  applying  this  correction  since  it  assumes 
horizontal  turbulence  stratification.  It  is  not  vEdid  when  large  contributions  of  dome 
or  mirror  induced  turbulence  dominate  the  situation. 

c.  Short  Exposure  Correction 

The  second  step  in  correcting  the  optical  data  is  concerned  with 
the  effects  of  exposure  time.  Fried  discusses  this  in  detail  [Ref.  17].  J.  Oldenettal 
[Ref.  18]  performed  additional  study  for  the  1.6  m  device.  Both  the  1.2m  and 
1 .6m  devices  measure  a  short  exposure  FWHM.  However,  inconsistencies  in  the 
1 .2m  Hughes  seeing  monitor  data,  made  comparison  impossible. 

The  acoustic  sounder  obtains  a  Cj*  profile  of  all  turbulence  it  detects 
and  provides  a  means  for  determining  a  corresponding  r^  value.  As  mentioned 
above,  the  1.6m  CIS  device  measures  this  same  turbulence  in  a  short  exposure 
manner.  In  essence,  this  means  that  the  imaging  system  eliminates  a  portion  of 
the  turblence  by  compensating  for  long  term  phase  tilts  in  the  light  received  from 
a  stellar  object.  Even  though  this  tip/tilt  correction  is  a  necessary  part  of  the  1 .6m 
compensated  imaging  system,  it  causes  the  device  to  detect  a  lower  level  of 


42 


turbulence  than  the  sounder.  This  lost,  long  exposure  turbulence  effect  must  be 
recovered  before  making  a  quantitative  comparison.  Table  1  contains  the  data 
taken  from  Reference  18  used  to  do  this.  The  second,  third,  and  fourth  entries 


TABLE  1.  CALCULATED  IMAGE  FWHM  FOR  LONG  AND  SHORT  EXPOSURES 
FOR  THE  1.6m  CIS  TELESCOPE.  [Ref.  18] 


To 

long  exp.  @  500nm 

short  exp.@  500nm 

(cm.) 

(arcseconds) 

(arcseconds) 

6 

1.69 

1.27 

10 

1.02 

.702 

14 

.734 

.446 

18 

.576 

.328 

22 

.472 

.234 

26 

.402 

.186 

30 

.352 

.142 

in  this  table  span  the  range  of  r^  values  measured.  A  least-squares  fit  of  the  1 .6 
m  long  exposure  r„  values  at  500  nm  to  the  reciprocal  of  the  short  exposure 
FWHM  values  gave 
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r^(cm)  =  2.975  + 


(10) 


4.925 
05 

with  a  correlation  coefficient  of  r  =  0.9999,  where  0s  is  the  short  exposure  FWHM 
corrected  for  zenith  angle  and  k  =  SOOnm.  Figures  23  through  25  contain  the 
final  results  of  these  conversion  procedures. 

Note  that  the  optical  r^  values  tend  to  lie  around  8  cm,  which  is 
considerably  below  the  12-15  cm  acoustic  results.  In  addition  the  optical  r^  values 
around  14:00  on  2  July  show  an  artificial  enhancement  produced  by  the  zenith 
angle  correction.  Both  of  these  results  suggest  that  there  was  considerable  dome 
or  mirror  turbulence  in  the  immediate  vicinity  of  the  1.6  m  (and  1.2  m)  devices. 
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Figure  23.  Acoustic  Sounder  Turbulence  Layering  Profile  Data  vs.  Corrected 
1.6m  CIS  Optical  Data  for  the  Evening  of  26-27  June,  1991. 
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Figure  24.  Acoustic  Sounder  Turbulence  Layering  Profile  Data  vs.  Corrected 
1.6m  CIS  Optical  Data  for  the  Evening  of  29-30  June,  1991. 
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Figure  25.  Acoustic  Sounder  Turbulence  Layering  Profile  Data  vs.  Corrected 
1.6m  CIS  Corrected  Optical  Data  fro  the  Evening  of  01-02  July,  1991. 
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4.  Wind  and  Atmospheric  Conditions 

As  discussed  in  the  Equipment  chapter,  one  of  the  recent  updates 
made  to  the  operating  software  was  a  Doppler  algorithm.  This  allows  for  the 
collection  of  horizontal  wind  speed  information  corresponding  to  the  turbulence 
under  investigation.  This  data  can  be  extremely  useful  in  localizing  regions  of 
high  turbulence  within  a  given  boundary  layer.  If  these  regions  are  low,  then  they 
can  be  completely  negated  through  the  use  of  an  elevated  platform.  This  is 
expected  to  be  the  case  for  smoothly  contoured  hilltop  locations.  However,  in  the 
case  of  the  AMOS  location,  the  mountaintop  is  not  smoothly  contoured.  In  fact, 
the  entire  facility  is  built  on  the  rim  of  an  ancient  volcano.  It  was  uncertain  what 
type  of  boundary  layer  would  be  formed  across  such  contours.  Doppler 
information  was  taken  during  the  Maui  visit  beginning  on  28  June  at  an 
inclination  angle  of  approximately  20  degrees  off  the  vertical  and  into  the 
direction  of  the  predominant  north-westerly  winds.  Figures  26,  27,  and  28  contain 
representative  examples  of  this  data.  The  altitude  cut-off  for  the  wind  data  is  the 
point  at  which  the  C/  signal  becomes  insufficient  to  be  considered  trustworthy. 
This  varies  constantly  and  is  a  direct  effect  of  the  wind  and  humidity  conditions. 
Appendix  C  contains  a  more  extensive  set  of  wind  profile  plots  collected  using 
the  Doppler  algorithm.  The  general  pattern  of  the  wind  speed  was  essentially  a 
constant  6  m/s  with  altitude,  although  a  10-20%  increase  appears  around  40-50 
m.  The  horizontal  wind  speeds  on  28-29  June  were  slightly  larger  being  around 
10  m/s. 
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In  addition  to  acoustic  echo-sounder  measurements,  a  VIZ  9000 
microsonde  system  made  measurements  of  humidity  and  temperature  profiles 
above  AMOS.  This  data  is  the  subject  of  a  Naval  Postgraduate  School  Master's 
thesis  being  written  by  C.  R.  Hoover  [Ref.  19].  Figures  29,  and  30  are 
representatives  of  the  low  altitude  portion  of  this  data  which  are  roughly 
coincident  in  time  with  the  wind  profiles  of  Figures  26  and  27. 


Figure  26.  Wind  Profile  for  Evening  of  28-29  June  1991. 
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Figure  28.  Wind  Profile  for  Evening  of  01-02  July  1991. 
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Figure  29.  Low-Altitude  Microsonde  Temperature  and  Humidity  Profiles  at 
AMOS  on  the  Evening  of  28-29  June  1991. 


Figure  30.  Low-Altitude  Microsonde  Temperature  and  Humidity  Profiles  at 
AMOS  on  the  Evening  of  29-30  June  1991. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  MISSION  OBJECTIVES 

1.  r,  vs.  Altitude 

Figures  19  through  22  imply  that  the  region  from  7.5-100  meters 
contains  a  significant  portion  of  the  optical  turbulence  encountered  above  AMOS. 
Additionally,  this  turbulence  occurs  in  stratified  layers.  This  suggests  that  the 
resolution  will  improve  by  positioning  the  telescope  on  an  elevated  tower.  Using 
Figures  23  and  24,  this  improvement  would  appear  to  be  on  the  order  of  30-40% 
for  a  15-25  m  tower.  This  value  was  found  by  comparing  the  existing  seeing  to 
that  available  at  15-25  m  and  neglecting  other  possible  turbulence  sources.  The 
conditions  on  the  evening  of  01  -02  July  were  anomalously  good  and  were  not 
considered  in  drawing  these  conclusions. 

Other  trends  noticed  in  the  sounder  results  are  that  the  seeing 
conditions  tend  to  improve  slowly  in  the  late  evening  hours  and  to  reach  a 
maximum  at  around  4  a.m.  each  morning.  This  suggests  that,  under  similar 
wind/weather  conditions,  a  window  of  optimum  viewing  may  be  predictable. 
Another  trend  is  the  rapid  decline  in  viewing  conditions  beginning  at 
approximately  5:30  a.m.  each  morning.  This  is  caused  by  convection  following 
sunrise.  The  interesting  point  is  that  it  seems  to  take  much  longer  for  the 
conditions  to  improve  after  sunset  than  it  does  for  them  to  decline  after  sunrise. 
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After  reviewing  the  entire  set  of  measured  wind  profiles,  nothing 
conclusive  can  be  said  about  the  wind  velocity  boundary  layer  above  AMOS. 
Characteristics  vary  considerably  from  plot  to  plot,  making  any  conclusions 
sketchy  at  best.  Since  there  is  no  apparent  drop  off  in  wind  speeds,  it  seems  that 
the  wind  shear  layer  extends  to  an  altitude  of  at  least  80  meters  above  AMOS, 
although  this  layer  is  only  optically  turbulent  in  its  lower  levels. 

2.  Site  Comparison 

As  demonstrated  in  Figures  1 7  and  1 8,  as  well  as  in  some  of  the  plots 
in  Appendix  B,  there  appears  to  be  no  significant  difference  between  the 
turbulence  profiles  at  Sites  A  and  G.  This  is  not  to  say  that  there  are  no 
differences,  but  that  this  data  set  does  not  reveal  any  detectable  differences. 
Surprisingly,  since  the  two  sites  were  at  different  heights  on  the  mountaintop,  this 
trend  suggests  that  the  turbulent  layer  follows  the  ground  contour  in  a  manner 
similar  to  observed  fog  advection.  However,  it  is  impossible  to  draw  any  definite 
quantitative  set  of  conclusions  based  only  on  this  data  set.  In  order  to  do  so 
would  require  not  only  additional  data,  but  also  data  gathered  during  different 
weather  and  wind  conditions.  Additionally,  data  from  tower  mounted  temperature 
probes  located  at  each  of  the  prospective  sites  would  be  useful  in  making  any 
such  conclusions. 
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B.  FUTURE  CONSIDERATIONS 

1.  Atmospheric  Turbulence 

Before  making  any  definite  conclusions  regarding  the  atmospheric 
turbulence  profile  above  AMOS,  additional  data  must  be  taken.  The  data 
collected  in  this  project  spanned  only  a  period  of  one  week  during  the  summer 
season.  To  conduct  an  accurate  survey,  it  is  recommended  that  additional 
measurements  be  taken  during  different  seasons  using  a  variety  of  equipment. 
This  equipment  should  include  not  only  an  acoustic  sounder,  but  balloon  sonde, 
temperature  probe,  and  optical  devices  as  well. 

2.  Artificial  Turbulence 

The  term  ‘artificial  turbulence’  refers  to  that  turbulence  caused  by 
devices  associated  with  the  telescope  being  used  for  observation.  These  include 
all  ventilation,  electrical,  and  building-related  heat  or  wind  turbulence  effects. 
Reference  12  presents  a  historical  study  of  such  effects  at  the  CFHT.  Figures  23 
and  24  imply  that  25-50%  of  the  optical  turbulence  could  be  dome  or  telescope 
induced.  The  instrumental  effects  of  the  1 .6m  telescope  and  the  CIS  were  not 
removed  from  the  optical  data  since  they  were  unknown.  If  the  optical  quality 
factors  were  available  and  included,  they  would  tend  to  improve  the  optical 
results  and  to  reduce  the  difference  with  the  sounder  data.  It  is  unlikely  that  this 
would  compensate  for  all  the  differences.  Measurements  of  dome-air  temperature 
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differences  and  turbulence  within  the  dome  and  over  the  mirrors  are  needed  to 
resolve  this  issue. 
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APPENDIX  A.  ACOUSTIC  ECHO-SOUNDER  OPERATING  CODE 

The  following  computer  code  is  used  to  operate  the  acoustic  echo-sounder 
discussed  in  this  paper.  It  is  written  in  Rocky  Mountain  Basic  and  is  for  use  on 
Hewlett-Packard  computers. 
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4  DEC  91,  PREUI0US:27  JUN  1991,  DLW 

The  computer  program  receives  information  from  an 'acoustic 
array  through  an  A-D  converter.  This  information  is  the 
returned  signal  of  an  acoustic  pulse  as  it  passes  through 
the  atmosphere.  The  data  is  then  used  to  display  the  return 
intensity  with  distance  as  a  function  of  time.  This  is  a 
short  range  device  (from  100m  to  150m). 


LIST  OF  UARIABLES 

A  -  The  area  of  the  receiver  (array) 

Amp$  -  Amplitude  for  HP33)4A 

Ap$ ,Frq$ ,Nm$ ,En$ ,Uo$  ,Hz$  -  strings  needed  to  set  the  HP3314A 
Counts  &  Cnt  -  the  string  and  real  representation  of 
Ct(«)  -  the  atmospheric  temperature  structure  parameter 

C0kelvin  -  the  speed  of  sound  at  0  degrees  celcius  in  m/s 
Ct_pts  -  the  number  of  traces  used  in  computing  Ct(*) 

Dat(*)  -  array  of  a-d  converter  output  after  sampling 
Disc_addre5s$  -  storage  location  of  data  file 
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End_pt  -  the  last  point  of  the  plot  vertically  taking 
Er  -  efficiency  of  coversion  of  acoustical  power 
to  electrical  power  on  the  recieve  side 

Et  -  the  efficiency  of  conversion  of  electrical 

power  to  acoustical  power  on  the  transmit  side 
range  correction  from  being  applied  to  data  very  near 
FreqS  &  Freq  -  the  frequency  input  of  the  HP3314A  in  kilohertz 
G  -  The  effective  aperture  factor 
Gain  -  the  electronic  gain  of  the  equipment 
Files  -  file  name  used  to  store  data 
Mrs  -  the  integer  hours 

I  &  J  &  K  -  counters  for  loops 

Ilvl  -  the  user  input  intensity  level  divisor;  this  value 
sets  the  screen  contrast  in  data  return 
Inc  -  The  vertical  increment  along  the  plot 
Kmax  -  total  number  of  block  averages  computed 
K3  -  the  wave  number  to  the  I /3rd  power 
M  -  the  counter  used  in  the  plot  label  routines 
Maxrng  -  the  maximum  range  of  the  echosounder  found 
Min  -  the  number  of  minutes 

Noise  -  a  running  total  of  the  noise  accumulated  at  the 
maximum  range;  used  to  find  and  remove  the 
average  noise 

Nrec  -  record  number  for  storage 

Ns  -  the  final  average  of  the  block  of  data  points;  also 
used  in  the  computation  of  the  noise  figure 
N$  Br  Num  string  and  real  representation  of  the  number  of 
cycles  in  the  transmitted  burst 
Numi  -  counter  for  computation  of  block  averages 

Num2  -  counter  for  the  number  of  block  averages  made 
Offset  -  the  computed  D-C  offset  for  the  system  prior 
Rec_num  -  real  and  integer  representation  of  the 

record  number  for  storage 

Plotnum  -  counter  used  to  insert  form  feeds  between  plots 
Pls_lng  -  the  pulse  length  of  the  burst 

Pnt  -  the  horizontal  position  on  the  Ct  plot 
Point<»)  -  real  number  representing  the  noise  and  range 
corrected  average  over  given  number  of  points 
Pt  -  The  computed  transmitted  power  to  the  array 

Otrhr  &  Qtrmin  -  keeps  track  of  the  passage  of  each 

15  minute  interval 

R  -  the  range  of  a  block  of  data  samples 
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810  I  Rdist  -  the  distance  traveled  in  one  tine  increment 

820  I 

830  •  by  multiplying  the  distance  per  time  increment 

840  !  by  the  number  of  time  increments 

850  !  Rspd  -  the  relative  speed  of  sound  that  the  echosounder 

880  I  sees  which  is  half  the  computed  speed 

870  I  Run_ave  -  the  running  sum  of  the  block  samples 

880  I 

890  !  Samfreq  -  the  user  input  sampling  frequency  desired 

900  !  Save_plt  &  Cntrl  -  on/off  toggles  used  to  determine  whether 

910  I  a  particular  run  will  be  saved  to  a  disc 

920  !  Sites  -  the  name  of  the  appropriate  site  of  data  collection 

930  !  Spd  -  the  speed  of  sound  in  air  computed  for  the  input 

940  !  temperature 

950  I  Temp  -  temperature  of  the  surroundings  in  degrees  celcius 

960  !  Timcnt  -  time  interval  between  data  samples  in 

970  I  Time  -  the  horizontal  position  of  the  trace  on  the  plot 

980  I  Timedist  -  the  horizontal  width  of  the  trace  on  the  plot 

990  !  horizontal  axis 

1000  I  to  data  acquisition 

1010  !  Times  -  the  string  required  by  the  A-D  converter  to 

1020  I  sample  at  the  desired  rate 

1030  !  Tminc  -  time  interval  between  A-D  samples 

1040  !  Trig  -  HP3314A  trigger  flag 

1050  !  Ypl  -  the  vertical  position  on  the  Ct  plot 

1080  •  Z  -  The  final  reduced  data  points  which  are  output  on  the  plot 

1070  !  Zeff  -  Effective  Impedence  for  input  loading 

1080  I  Zimp  -  the  speaker  array  impedence 

1090  I  Zones  -  the  appropriate  time  zone  the  operator  desires 
1100  I 

1110  OPTION  BASE  1 
1120  I 

1130  I  initialize  the  arrays  &  set  dimensions 
1140  I  declare  all  integer  variables 
1150  ! 

1180  DIM  Disc_address$[20]  .File1St30]  ,File2S[30] ,D1<500  )  ,Ct(250  >  .Point (300 ) 
1170  DIM  Doppler<  300 ) .Vdopl 300  >  ,V_num( 300  > ,Davg( 2000 ) ,Ct_dop( 600  ) 

1180  INTESER  Hr  ,I  .Inurn  ,J  ,K  ,Kmax  ,K1  ,K2  ,M,Ncnt ,Nbin ,Num ,Num1  ,Nun2 
1190  INTEGER  Plotnum ,Pr int_key 

1200  DIM  Dat< 18200)  !  FOR  20  KHZ 

1210  ! 

1220  I  initialization  rout ines ....  set  time,  set  HP3314A  Function 
1230  !  Generator 

1240  ! 

1250  Rstrt:CALL  Freq_ini t < FreqS  ,NS  ,V_out  ) 

1260  CALL  Init_ad200 

1270  CALL  Set  time(ZoneS) 
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INPUT  "SITE  NAME  " .SiteS 


keyboard  set_iip 


sets  labels  on  the  computer  function 
keys 


OUTPUT  KBD? “SCRATCH  KEYE* 
CONTROL  2.2 i1 


!  CLEAR  THE  KEYS 


1350 

ON 

KEY 

0 

LABEL 

"NEW 

TEMP"  GOTO  Speed 

1360 

ON 

KEY 

I 

LABEL 

"FILTER 

GAIN"  GOTO  Set_gain 

1370 

ON 

KEY 

2 

LABEL 

"PRINT 

TRACE"  GOTO  Plt_dmp 

1380 

ON 

KEY 

3 

LABEL 

“COLOR 

SCALE"  GOTO  Doppler_f5 

1390 

ON 

KEY 

4 

LABEL 

"PLOT 

RANGE"  GOTO  Plot_range 

1400 

ON 

KEY 

5 

LABEL 

"RESTART" 

GOTO  Rstrt 

1410 

ON 

KEY 

9 

LABEL 

"QUIT"  GOTO  Quit 

GINIT 


I  Set  constants 


Area= . 0866 

Blue=2/3 

Count$=" 18200” 

C0kelvin=331 .6 

Di sc_addres5$= " : ,70 

Er=.5 

Et=.5 

Filter_gain=10 
Ful l_scale*5 
6=.4 

Green= 1 /3 

Ilvl=.01 

Maxsec=900 

Max_rec I =96 

Max_recZ=22000 

Plotnum=1 

Pltrng=l 00 

Power_amp_gain=20 . 6 

Rin=100 

Samf req=20000 

Samave=0 

I Z  imp* 12.12 

Uloss=2 . 2 

Zoff-12.94 

Zimp=1 9.5 


!  Array  area  m^Z  0 . 00456/dri ver 
I  Color  value  for  pure  blue 
I  Totalf  A-D  samples  11200  or  18200 
I  Speed  of  sound  at  0  C 


704 .0  .0' 


!  Receive  efficiency 
!  Transmit  efficiency 

I  Uoltage  gain  for  filter  gain  of  20  db 
!  Doppler  full  scale  (red  or  blue) 

I  average  antenna  gain  factor 
!  Ualue  for  green 
!  CT  threshold  for  full  scale 

!  number  of  seconds  in  a  15  min  plot 
'  96  files  in  24  hrs ,  1  every  15  min 
!  tCt2  traces  stored  in  24  hours  @  45ec 

!  Max  range  on  plot 

•  Power  amp  voltage  gain 

!  preamp  input  impedance ,  ohms 

•  sample  frequency 

•5x5  square  array  impedance 
I  Voltage  loss  in  line  and  TR  switch 
i  effective  Z  for  100  ohm  input 
I  19  ELEMENT  HEXAGONAL  ARRAY  Z 


I  Derived  constants 

INPUT  "ENTER  PRE-AMP  GAIN:  987.1906  or  3895" .Preamp_gain 
INPUT  "SET  FILTER  GAIN  ON  PREAMP  AND  ENTER  SETTING  ( 0 . 1 0 ,20  ,30  )“ .Db 
PRINT  "FILTER  GAIN  =  " .Db 
F 1 1 ter_gain= 1 0" ( Db/20  > 
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1760 
1770 
1780 
1790 
1800 
1810 
1820 
1830 
1840 
1850 
1860 
1870 
1880 
1890 
1900 
1910 
1920 
1  950 
1  940 
1  950 
1960 
1  970 
1  980 
1990 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2080 
2090 
2100 
2110 
2120 
2130 
2140 
2150 
2160 
2170 
2180 
2190 
2200 
2210 
2220 
2230 
2240 


I 

'Jtpk=i^'_out  ♦Pouier_amp_gain 

Pt  =Et  *  ( t'tpk-Uloss  >  "2/(  2 .  *2  iMp  )  I  RMS  acoustical  power  output 

Gain=Preamp_gain*F i 1 ter_gain 

I 

I 


Sfrq;  I 

6CLEAR 


I 


I 


Tpn:nc=l  .  /Sanfreq 

Tincnt=Tninc/ 1 .E-9 

Time$="TIME  (  Timcnt  ) 

Cnt='JAL(  Counts  ) 

Navg_pt5*2000 

Ncnt'Cnt 

Nun-UAKNS  / 

Freq=0AL<  FrcqS )* 1 000 
Fu 1 l_scale=S 
Tau=Num/Freq 
Nbin=Safnfreq*Tau 
K5ave=INT(Ncnt/Nbin) 
Navg=Ncnt -Navg_p  t  s+ 1 
Nb 1 n2=Nb i n/2 
Numi =Nbin-1 
Num2=Ncnt-Nbin+1 
Nun_pts=Nbin+Nbin  MOD  4 


Delta  t  in  nsec 

Total  *  of  A-D  samples 
t  points  used  to  compute  offset 

t  cycles  transmitted  in  burst 
Transmitted  frequency 
Doppler  full  red  or  blue  shift 
Pulse  length  ( sec  ) 

#  of  counts  per  range  bin 

#  of  range  bins 


REDIM  Ct ( Ksave ) ,Poi nt ( Ksave  )  ,Dopp ler( Ksave ) ,D 1 ( 2»K5ave+ 1 0  )  ,Udop<  Ksave  ) 
REDIM  l'_nun''  Ksave  )  ,Ct_dop<  2»K5ave  > 


CALI  Fi  le_ini  t  (  Di  sc_address$  ,Max_rec  1  ,Ma.v_rec2  .Mlast  1  .NlastC  , Ksave  ,Fi  1 
ASSIGN  eFilel  TO  FilelS 
ASSIGN  @File2  TO  File2$ 

Nrec_opt=Nla5t2 


computation  of  the  speed  of  sound  at  a  given  temp  and  the 
range  of  detection  of  the  device 


Speed: INPUT  “Enter  the  temperature  (Celsius)  “.Temp 
INPUT  “Enter  the  relative  humidity  in  /(“  ,Rel 
INPUT  "Enter  the  atmospheric  pressure  in  mb"  .Atmo5_pres 
Spd“C0kel vin*( SQR( 1 .  +  ( Temp/273.  > )  > 

R5pd“Spd/2 . 

Rdi5t=Tminc*R5pd 

Maxrng=Cnt  *Rdi st 

Lambda*Spd/Freq 

K3*< (2.*PI )/Lambda  )"( 1 . /3.  ) 

r l5_lng“(R5pd*Num  )/Freq 
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2250  I 

2260  CALL  At tenuat ion( Freq  .Temp  ,Re 1 ,Atno5_pres  ,At ten > 

2270  At ten2“2 . *At ten 

2280  ! 


2290  Start_plot :  I 
2300  ! 

2310  !  Computation  of  the  D-C  offset 

2320  ! 

2330  Noi5e“0 
2340  Off5et=0 


2350 

2360 

2370 

2380 

2390 

2400 

2410 

2420 

2430 

2440 

2450 

2460 

2470 

2480 

2490 

2500 

2510 

2520 

2530 

2540 

2550 

2560 

2570 

2580 

2590 

2600 

2610 

2620 

2630 

2640 

2650 

2660 

2670 

2680 

2690 

2700 

2710 

2720 

2730 


Of  f set_f lag=0 
Odoppler=Spd/< 2*PI  ) 

Color_5cale*1  /(  3*Ful l_5cale  ) 

Uol 1 5_per_bi t  =5./<2048.«6ain*Rin/(Rin+Zeff  >  > 

!  Compute  CT  calibration  constant 
Ct_cal=Uol  t  s_per_b  it'‘2/2imp 
Ct_cal=Ct_cal*(  Temp +273  )''2/(  .0039*K3) 

Ct_cal=Ct_cal/<  Area*6*Pt*Er*Pi5_lng ) 

I 

Trig*0 

CALL  Read_ad200<  Oat ( *  )  ,T ime$  .Trig  ) 

Trig=1 

I 

PRINT  "COMPUTING  D-C  OFFSET  AND  NOISE  BACKGROUND" 

Of f set=SUM( Dat  )/Cnt 
MAT  Dat*  Dat-<Offset) 

I 

I  Compute  real  and  imaginary  noise  levels 

I 

K=0 

Zreal=0. 

Z inag=0 . 

FOR  1=1  TO  Num2  STEP  Nbin 
K=K  +  1 

Inum«I+Num1 

FOR  J*I  TO  Inum  STEP  4 
A=Dat  <  J  ) 

B*Dat  <  J+ 1  ) 

C-Dat<  J+2  ) 

D*Dat ( J+3  ) 

Zreal=A*C+B*D+ZreaI 
Zimag=-A*D+B*C+Zimag 
NEXT  J 
NEXT  I 

Zrnci 5e=Zreal /K 
Zi*Zimag/K 

Z 1  no  1  se*Zrno 1 5e*TAN( 1 /Udopp ler  )  '  +1m/5  noise  offset 

PRINT  Zrnoi se , Z 1  .Zinoise 
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2740  !  Compute  total  noise 

2750  ! 

2760  MAT  Dat*  Oat  .  Dat 

2770  Noise=SUM(Dat )/Cnt 

2780  ! 

2790  I  set  up  the  plot 
2800  ! 

2810  FOR  Nrec=Nlast1  TO  Max_rec1 

2820  CALL  Plot_setup(  Maxsec ,5ite$ .Maxrng , PI trng , Zones  ) 

2830  MAT  Ct“  <0. )  !  Zero  arrays 

2840  MAT  Ooppler*  <0.  ) 

2850  MAT  Wdop-  <0. ) 

2860  MAT  U_num=  (0.  ) 

2870  Ct_pt5=0. 

2880  Nsec_la5t=0 

2890  OUTPUT  KB0;"L‘'; 

2900  ! 

2910  >  Compute  15  minute  noise  level 

2920  ! 

2930  Trig=0 

2940  CALL  Read_ad200( Oat (♦)  .Times .Trig  ) 

2950  Trig=I 

2960  ! 

2970  Of fset  =  SUM(Oat  )/Cnt 

2980  MAT  Dat*  Dat  .  Dat 

2990  Neu_noise=SUM(  Oat >/Cnt-0f f set*0f f set  I  Noise  variance 

3000  IF  Noise=0  THEN 

3010  Noise=Neu_noise 

3020  ELSE 

3030  IF  New_noi se^Z . *Noi se  THEN  Noise*. 5*Noise+.5*Neu_noi5e 

3040  END  IF 

3050  Signal_thresh*2 . ‘Noise  I  Signal  threshold  for  winds 

3060  I 

3070  Read_sig:  '  read  the  A-D  converter 

3080  I 

3090  I  find  the  t  seconds  in  the  fifteen  minute  block 
3100  ! 

3110  T0-TIMEDATE 

3120  Tt=(T0  MOO  86400)  i  $  sec  past  midnight 

3130  N5ec  =  (T1  MOD  3600  MOD  Maxsec  )  i  #  sec  within  15  min  block 

3140  Delta_t ime=Nsec-Nsec_last  !  Time  interval  between  pulses 

3150  Nsec_la5t*Nsec  '  Old  #  seconds 

3160  IF  Delta_time<0  THEN  SOTO  Plt_dmp  '  Print  plot  after  15  minutes 

3170  Del ta_t ime*Del ta_t ime/2+2 

3180  I 

3190  CALL  Read_ad200( Dat (*>  .Times  .Trig  ) 

3200  I 

3210  MAT  Oavg*  Dat < Navg : Ncnt  ) 


62 


3220 

3230 

3240 

3250 

32G0 

3270 

3280 

3290 

3300 

3310 

3320 

3330 

3340 

3350 

3360 

3370 

3380 

3390 

3400 

3410 

3420 

3430 

3440 

3450 

34B0 

3470 

3480 

3490 

3500 

3510 

3520 

3530 

3540 

3550 

3580 

3570 

3580 

3590 

3600 

3610 

3620 

3630 

3640 

3650 

3660 

3670 

3680 

3690 


Of f set=SUM(  Davg )/NavQ_pts 
MAT  Oat*  Dat-(Off5et) 

MAT  Odop*  (0.  ) 

! 

K=0 

ALPHA  OFF 

FOR  I“1  TO  Nut*i2  STEP  Nbin 

R=( I+Nbin2  )*Rdi5t  !  Mean  range  of  range  bin 

K=K+1 

Rau)_5ig=0. 

Inun=I+Nun1 

2real=0. 

Zirnag=0. 

FOR  J=I  TO  Inun  STEP  4 
A=Dat( J  ) 

B=0at( J  +  1  ) 

C-Dat<  J+2  ) 

D=Dat( J+3  ) 

Zreal=A«C+B*0+2real 
Z inag=-A*D+B*C+Zinag 
Rau_5ig=Raw_5 ig+A*A+B*8+C*C+D*D 
NEXT  J~ 

Raw_5iQ=Raw_sig/Nun_pt5 

Poi nt  <  K  )=Ct_caI *R»R*( Raw_5ig-Noise  >*EXP( R*At  ten2  > 

IF  Raw_sig>Signal_thresh  THEN 

Udop<  K  )*ATN<  <  Zinag-Zinoise  >/<  Zreal-Zr noise  )  )*Udoppler 
V_nun(  K  )*y_nu(n(  K  H I 
END~IF 
NEXT  I 

Ki*ia,x=K 

I 

I 

'  Add  data  to  Ct  average  if  the  noise  is  sufficiently  low 
I 

IF  Raw_sig<20. *Noise  THEN  I  Add  Points  TO  Ct 

MAT  Ct=  Ct+Point 
MAT  Ooppler=  Doppler+Wdop 
Ct_pt5*Ct_pt  5+ 1 
ELSE 

OISP  “HIGH  BACKGROUND  NOISE-  SIGNAL  REJECTED- 
END  IF 

I 

'  Plotting  of  the  data 

I 

WINDOU  0  .Maxsec  ,0 ,P 1 trng 
GRAPHICS  ON 

Inc*Nbin*Rdist  '  Uertical  increment 
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3700  IF  Delta_t ine> 1 0  THEN  Delta_ti«e“10  !  account  for  delays 
3710  Color_5cale= 1 /( 3*Ful l_5cale  )  !  Set  full  scale  color 

3720  ! 

3730  !  compute  the  intensity  of  the  return,  move  to  the  proper 
3740  I  coordinates  and  plot  the  appropriate  dithered  block 
37B0  ! 

3760  FOR  K*1  TO  Kmax 
3770  Pointk=Point<K ) 

3780  IF  Pointk>0.  THEN 

3790  2*SQR(Pointk/Ilvl > 

3800  IF  Z>1 .  THEN  Z=1  . 

3810  ELSE 

3820  Z=0 

3830  END  IF 

3840  0color=0dop(  K )*Color_scale+6reen 

3850  IF  Ocolor<0  THEN  Ucolor=0  I  Red 

3860  IF  Ucolor>Blue  THEN  Ucolor=Blue  '  Blue 

3B70  AREA  COLOR  Ocolor  .  1  .Z 

3880  Range=<K-.5  )*Inc 

3890  MOVE  Nsec  .Range 

3900  RECTANGLE  Delta_t ime  .Inc  .FILL 

3910  NEXT  K 

3920  I 

3930  Point ( 1 )=T0 

3940  MAT  Ct_dop(  1  :Ksave  >=•  Point 

3950  MAT  Ct_dop< Ksave+I :2*K5ave >=  Wdop 

3960  OUTPUT  @Fi le2 ,Nrec_opt i Ct_dop( •  ) 

3970  Nrec_opt=Nrec_opt+ 1 
3980  I 

3990  GOTO  Read_5ig 
4000  I 

4010  !  graphics  dump  of  plot  after  15  minute  intervals 

4020  I 

4030  Plt_dmp:  PRINTER  IS  701 

4040  PRINT  ■■ 

4050  PRINT 

4060  PRINT  " 

4070  Plotnum=Plotnum+1 

4080  DUMP  GRAPHICS  #701 

4090  I 

4100  MAT  Ct=  Ct/<Ct_pts  ) 

4110  FOR  I* 1  TO  Kmax 

4120  IF  Ct(I)<1,E-6  THEN  Ct(I)=1.E-6 

4130  IF  U_num(I)>0.  THEN  Doppler( I >=Doppler( I  )/U_num( I  ) 

4140  NEXT  I 

4150  I 

4160  CALL  Ct_plot ( Ct < *  )  ,P1 1 rng  .Maxrng  .Kmax ,Si te$  ) 

4170  I  Save  the  Ct  data 
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!  ZERO  ARRAY 
I  T0  =  TIMEDATE 


4180  MAT  01-  (0.  > 

4190  D1<1>=T0 

4200  01  <  2  >=Atroo5_pres 

4210  D1<3>=Tenp 

4220  D1<4>=Rel 

4230  01(S)=Atten 

4240  D1(6)-Inc 

4250  D1<7)-Ct_cal 

4260  D1<8)-Gain 

4270  D1(9>=Kfnax 

4280  Dl(10>=Off5et 

4290  ! 

4300  Ki=10+K5ave 

4310  FOR  1=1  TO  Kmax 

4320  01 ( 1  +  10  )=Ct( I  ) 

4330  01 ( I+Ki  )=Doppler( I  ) 

4340  NEXT  I 

4350  ! 

4360  OUTPUT  @F:le1  .Nrec;D1(*) 

4370  I 

4380  I  start  the  next  15  minute  plot 

4390  I 

4400  next  Nrec 
4410  ! 

4420  Doppler_fs:  ' 

4430  INPUT  "Enter  Doppler  full  scale  in  m/s.  (5  is  normal )" ,Ful l_5cale 
4440  GOTO  Read_5ig 

4450  Plot_range:  ' 

4460  INPUT  "Enter  CT  plot  full  scale  10t  or  150  m” ,Plot_rng 
4470  PI trng=INT( Plot_rng/50  )*50 

4480  GOTO  Start_plot 

4490  Set_gain:  I 

4500  INPUT  "Enter  filter  gain,  0,  10.  20,  30,  40  OB"  ,Db 

4510  IF  Db=0  OR  Db=10  OR  Db=20  OR  Db-30  OR  Db=40  THEN 

4520  Fi 1 ter_ga in= 1 0" ( Db/20  ) 

4530  ELSE 

4540  PRINT  "URONG  GAIN  SETTING” 

4550  GOTO  Set_gain 

4560  END  IF 

4570  Gain=Fi lter_gain«Preamp_9Bin 

4580  GOTO  Start_plot 

4590  Ou  1 1 :  * 

4600  PRINT  "Output  has  been  written  to  "iFileS 

4610  ASSIGN  @Fiiel  TO  ♦ 

4620  END 

4630  I 

4540  I 

4650  I  SUBROUTINE  SECTION 
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46E0  • 

4670  SUB  Freq_ini t( FreqS ,N$  ,U_out  ) 

4680  I 
4690  I 

4700  I  SETUP  OF  THE  HP3314A  FUNCTION  GENERATOR 
4710  I 
4720  ! 

4730  Ap$="AP" 

4740  Frq$="FR" 

4750  Ni*i$-"NM" 

4760  En$="EN" 

4770  Uo$=“UO‘' 

4780  H2$=“KZ" 

4790  INPUT  "FREQUENCY  DESIRED  (kHZ)  (5  RECOMMENDED )"  .Freq$ 

4800  INPUT  "AMPLITUDE  DESIRED  (U  .  .2.0U  MAX)  "  .Anp$ 

4810  IF  UAL< Anp$  )>2.0  THEN 

4820  Apnp$*"2.0" 

4830  PRINT  “AMPLITUDE  OF  FUNCTION  GENERATOR  IS  2.0  U" 

4840  END  IF 

4850  U_out  =  UAL( Amps  ) 

4860  INPUT  "NUMBER  OF  CYCLES  PER  BURST  (INTEGER)  <100  RECOMMENDED  )" ,N$ 

4870  OUTPUT  707i"MO3" 

4880  OUTPUT  707; "SR2" 

4890  OUTPUT  707;Ap$&Arnp*&Uo$&Frq$&Freq$&H2$&Nm$&N*&En$ 

4900  SUBEND 
4910  ! 

4920  I 

4930  SUB  Init_ad200 
4940  I 
4950  I 

4960  !  INITIALIZATION  OF  THE  A-D  CONVERTER 

4970  ! 

4980  ! 

4990  Ad_5el_code= 1 7 

5000  Dumny“READIO( Ad_5el_code ,3  ) 

5010  WRITEIO  Ad_5el_code  ,0;0 

5020  CONTROL  Ad_5el_code  ,0; I 

5030  SUBEND 
5040  ! 

5050  I 

5060  SUB  Set_t  ii*ie<  Zones  ) 

5070  ! 

5080  I 

5090  !  SET  THE  TIME  DATE  RECORDER 

5100  ! 

5110  ! 

5120  PRINT  "WHAT  TIME  REFERENCE  ARE  YOU  USING*?  INPUT:" 

5130  PRINT  “  1  FOR  UNIVERSAL  TIME" 
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5140 

5150 

5160 

5170 

5180 

5190 

5200 

5210 

5220 

5230 

5240 

5250 

5260 

5270 


PRINT  "  2  FOR  LOCAL  TIME" 

PRINT  ■■  3  FOR  YOUR  OWN  CLASSIFICATION- 

INPUT  K 
IF  K-2  THEN 

Zone$= "( LOCAL  )" 

ELSE 

IF  K=3  THEN 

INPUT  "WHAT  IS  YOUR  TIME  REFERENCE" ,Zone$ 

ELSE 

Zone*=*"(UTC  )" 

END  IF 
END  IF 

INPUT  "Do  you  want  to  reset  the  clock  7”  OS 
IF  Q$=“YES"  OR  Q$="Y“  THEN 


5280 

5290 

5300 

5310 

5320 

5330 

5340 

5350 

5360 

5370 

5380 

5390 

5400 

5410 

5420 

5430 

5440 

5450 

5460 

5470 

5480 

5490 

5500 

5510 

5520 

5530 

5540 

5550 

5560 

5570 

5580 

5590 

5600 

5610 


INPUT  “ENTER  “"DD  MMM  YYYY . ,Date$ 

INPUT  “ENTER  ““HR:MIN:SC . .TineS 

SET  TIMEDATE  DATE( OateS >+TlME( TineS > 
PRINT  DATE$( TIMEDATE  )  .TiMESt  TIMEDATE  ) 
END  IF 

SUBEND 


SUB  Read_ad200( Dat  <  *  )  .TimeS  .Trig  ) 
INTEGER  D1( 1:3200)  BUFFER 
INTEGER  D2( 1:15000)  BUFFER 
Count  1 S* "3200“ 

Count2$=" 15000" 


INFOTEK  A-D  ROUTINE  SET  UP  FOR  EXTERNAL  TRIGGER 


Ad_sel_code= 1 7 

I 

!  INITIALIZATION  OF  THE  A-D  CONVERTER 
( 

OUTPUT  Ad_5el_code! "RESET" 

OUTPUT  Ad_sel_code; "INTERNAL" ."COUNT  "&Count 1  $ . "HOLDON" 
OUTPUT  Ad_sel_code; "DELAYON" ."SELECT  Isl  end"  .TimeS 
OUTPUT  Ad_5el_code; "STATUS" 

ENTER  Ad_5el_code;Re5p$ 

IF  RespS«=" - "  THEN 

ASSIGN  @Ad200  TO  Ad_sel  code;W0RD 

I 

'  triggering  of  the  HP3341A 

I 

ASSIGN  @Buf1  TO  BUFFER  D1(*) 

ASSIGN  @Buf2  TO  BUFFER  D2(*) 
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5620 

5630 

5640 

5650 

5660 

5670 

5680 

5690 

5700 

5710 

5720 

5730 

5740 

5750 

5760 

5770 

5780 

5790 

5800 

5810 

5820 

5830 

5840 

5850 

5860 

5870 

5880 

5890 

5900 

5910 

5920 

5930 

5940 

5950 

5960 

5970 

5980 

5990 

6000 

6010 

6020 

6030 

6040 

6050 

6060 

6070 

6080 

6090 


IF  Trig=1  THEN  TRIGGER  707 
TRANSFER  @Ad200  TO  §Buf  1  ;IJAIT 

OUTPUT  Ad_5el_code; "COUNT  "&Count2S  .  "SELECT  IslO  end’ 
TRANSFER  @Ad200  TO  @Buf2;WAIT 
MAT  Dat<3201 : 1 8200  )*  D2 
MAT  Dat»  Dat*< . 1 ) 

MAT  Dat< 1 :3200  )=  01 

ELSE 

PRINT  "ERROR  DURING  INITIALIZATION  =";Re5p$ 

END  IF 
SUBEND 
! 

! 

SUB  Plot_5etup( Maxsec ,Site$  .Maxrng  .Pltrng  ,Zone$  ) 

I 

I 

I  SET-UP  OF  THE  TIME  PLOT  ON  THE  CRT 

I 

I 

IF  Pltrng=150  THEN 
Nun_label s=3 
Nun_t icks*! 5 
END  IF 

IF  Pltrng=100  THEN 
Nurn_label  5=4 
Nun_t ick5=20 
END  IF 

GCLEAR 

GRAPHICS  ON 

LINE  TYPE  1 

UIEWPORT  15,120,15.80 

WINDOW  0  , Maxsec  ,0 ,P 1 trng 

AXES  Max  sec/15  ,P1  trng/Num_t  icks  , 0,0, Max  sec /3  ,P1  trng/Nurn_label  s 
CLIP  OFF 
CSIZE  4, .6 
LORG  6 

Tl-TIMEDATE  MOD  86400 
Hrs=T1  DIU  3600 
T2-T1  MOD  3600 
Min=T2  DIU  60 
Otrhr=Min  DIU  15 

FOR  M=0  TO  Maxsec  STEP  INT( Maxsec/S > 

MOUE  M,-Pltrng/45 
Qtrnin*0trhr«15+<M*3/Max5ec  )*5 
IF  Otrmin=60  THEN 
Qtrmin=0 
Hrs-Hrs+I 
END  IF 
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6100  LABEL  USING  “DD  .A  ,ZZ” ;Hr5 ; " : " :Qtrmin 

6110  NEXT  M 

6120  MOUE  Maxsec/2,-15 

6130  LABEL  "TIME  “&Zone$ 

6140  I 

6150  I  LABEL  ORDINATE 
6160  I 

6170  LORG  8 

6180  FOR  M=0  TO  PltroQ  STEP  INT< Pltrng/Num.labels  ) 

6190  MOUE  0.M 

6200  LABEL  M 

6210  NEXT  M 

6220  LOIR  PI/2 

6230  LORG  6 

6240  MOUE  -Maxsec/7 .Pltrng/2 

6250  LABEL  "RANGE  IN  METERS" 

6260  I 

6270  I  TITLE 
6280  I 

6290  loir  0 

6300  LORG  4 

6310  MOUE  Max5ec/2 ,Pltrng+3 

6320  LABEL  OATE$< TIMEDATE ) ; "  ";Site$ 

6330  CLIP  ON 

6340  SUBEND 
6350  ! 

6360  ! 

6370  SUB  File_init(Disc_addres5$,Nrecl  ,Nrec2  .NIast 1 .Niast2 .Ksave  .Fi le 1  $ .Fi le2$ 
6380  ! 

6390  I 

6400  I  CREATE  THE  STORAGE  FILE  ON  THE  DISC  FOR 
6410  I  THE  REDUCED  DATA 

6420  I 
6430  I 

6440  INPUT  "ENTER  THE  REDUCED  DATA  OUTPUT  FILENAME  ",File1$ 

6450  INPUT  "Is  this  the  first  entry  m  the  file'i’",0$ 

6460  Fi  le2S=Fi  le  1 4 [  1  ,7  ]8t"EXT"&0i5C_addre55$ 

6470  P i le 1 $=Fi 1 e 1 $&Di sc_addre55$ 

6480  IF  0$="YES"  OR  Q$="Y"  THEN 

6490  CREATE  BOAT  Fi le1$  .Nrecl  ,8«( 2*K5ave+10  ) 

6500  CREATE  BOAT  Fi le2$  .Nrec2  ,8*2*K5ave 

6510  END  IF 

6520  ASSIGN  @File1  TO  FileU 

6530  STATUS  ©Filel ,7;Nla5t1  i  next  record  in  file  1 

6540  ASSIGN  @File2  TO  File2$ 

6550  STATUS  @Fi le2  ,7iNla5t2  !  next  record  in  file  2 

6560  I 
6570  SUBEND 
6580  I 
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G590  ' 

GB00  SUB  Ct_plot(Ct( *) .Pltrng  .Maxrng  .INTEGER  Knax.SiteS) 

BGI0  I  CT-'2  PLOT  ROUTINE 

GG20  •  Pitrng  is  the  naxipiun  range  in  50  rn  increments  (  1 50n  typically) 
BG30  I  Maxrng  is  the  actual  maximum  range  that  depends  on  temperature 
GG40  I  Kmax  -  the  maximum  number  of  Ct  dat  points  CtfKmax). 

6G50  6CLEAR 

GBB0  IF  Pltrng=I50  THEN 

GG70  Num_label5*3 

BG80  Num_ticks=15 

6G90  END  IF 

6700  IF  Pltrng=100  THEN 

6710  Num_labels=4 

6720  Num_tick5=20 

6730  END  IF 

6740  Decades=5 

6750  UIEWPORT  15,120.15,80 

6760  WINDOW  0 .Decades .0 ,P1 trng 

6770  AXES  1  ,P1 trng/Num_t ick s  .0  .0 , 1  ,P1 trng/Num_label s 

6780  CLIP  OFF 

6790  CSI2E  3. .9 

6800  LORG  6 

6810  FOR  M=0  TO  Decades 

6820  MOUE  M.-Pltrng/4S 

6830  LABEL  (M-Decades) 

6840  NEXT  M 

6850  MOUE  Decades/2,-15 

6860  CSI2E  4 

6870  LABEL  "LOG  OF  CT''2'‘ 

6880  LORG  8 

6890  FOR  M=0  TO  Pitrng  STEP  INT( Pltrng/Num_labels  ) 

6900  MOUE  0.M 

6910  LABEL  M 

6920  NEXT  M 

6930  LOIR  PI/2 

6940  LORG  6 

6950  MOUE  -Decades/T ,Pltrng/2 

6960  CSI2E  4 

6970  LABEL  "RANGE  IN  METERS” 

6980  LOIR  0 

6990  LORG  4 

7000  MOUE  Decades/2 , PI trng+3 

7010  LABEL  “TIME  AUERAGED  CJ-'Z  ”;Site$ 

7020  CLIP  ON 

7030  MOUE  Decades  ,0 

7040  Inc=Maxrng/Kmax 

7050  FOR  1=1  TO  Kmax 

70S0  Ypl-Inc*(I-.5 ) 
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7070  Pnt=LGT(Ct(I  )  )+Decades 

7080  IF  Pnt>Decades  THEN  Pni^Decades 

7090  IF  Pnt<0  THEN  Pnt=0 

7100  DRAW  Pnt  ,Ypl 

71 10  NEXT  I 

7120  DUMP  GRAPHICS  #701 

7130  PRINT  •' 

7140  PRINTER  IS  CRT 

7150  GCLEAR 

71S0  SUBEND 

7170  SUB  At tenuat ion< Freq  .Temp  ,Rel ,Atmo5_pres .Alpha  ) 

7180  I  11  June  1991:  DLW 
7190  I 

7200  I  This  subprogram  calculates  the  attenuation  coefficient 

7210  •  of  sound  in  air  based  upon  equations  in  H.  BASS, 

7220  !  Acoustical.  Soc .  Am.(xxx  1991) 

7230  I 

7240  !  Alpha  -  Intensiy  attenuation  coefficient  of  sound 

7250  •  (alpha  is  normally  calculated  for  pressure) 

7260  I  Atmos_pres  -  atmospheric  pressure  in  mb 
7270  I  Es  -  Saturation  vapor  pressure 
7280  •  Freq  -  sound  frequency  in  Hz 

7290  !  H  -  absolute  humidity  (mb) 

7300  I  Rel  -  relative  humidity  in  % 

7310  I  Temp  -  Celcius  temperature 

7320  I 

7330  !  Saturation  vapor  pressure  using  Hooper  formulation 

7340  T=Temp 

7350  A0» 1.3521 

7360  Al=1.6369E-2 

7370  A2-3.1794E-5 

7380  A3=-l .4892E-7 

7390  Es=(A0+T*(A1+T*(A2+T*A3)))"6 

7400  H=Rel*E5/AtmQ5_pres  !  Absolute  humidity  in  'L 

7410  I 

7420  I  H.  Bass  Attenuation  calculation 
7430  ) 

7440  F2=Freq*Freq 

7450  P50=Atmos_pre5/ 1 0 1 3 . 25  !  Relative  pressure 

7460  Tk*Temp+273. 15  '  KELVIN  TEMPERATURE 

7470  Tr-293.15  !  REFERENCE  TEMPERATURE  K 

7480  Tkr=Tk/Tr 
7490  I 

7500  •  Compute  nitrogen  and  oxygen  relaxation  frequencies 

7510  Frn=Pso/SOR(Tkr  )*( 9+280»H*EXP( -4. 179*(  (  1  /Tkr )'. 3333333- 1  > ) ) 

7520  Fro=P5O*(24+4.04E+4*H*( .02+H)/( .391+H  > ) 

7530  ! 

7540  Alpha-F2*( 1 .84E-1 1 ♦SQRI Tkr  )/Pso+Tkr " ( -2 . 5 )* ( 1 . 278E-2*EXP( -2239 . 1 /Tk  )/( Fro+ 
7550  Alpha=2*Alpha  I  convert  to  nepers/m  (intensity) 

7560  SUBEND 


71 


APPENDIX  B.  ACOUSTIC  SOUNDER  DATA  FROM  AMOS 

The  following  is  a  representative  set  of  acoustic  echo-sounder  profiles 
collected  at  AMOS  during  June- July,  1991.  (Note:  The  profiles  need  to  be 
multiplied  by  a  factor  of  3.8  to  be  calibrated.) 
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Figure  31.  Echo-sounder  Turbulence  Profile,  26  June,  1991. 
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Figure  32.  Echo-sounder  Turbulence  Profile,  26  June,  1991. 
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Figure  33.  Echo-sounder  Turbulence  Profile,  26  June,  1991. 


75 


RRNGE  IN  METERS  RRNGE  IN  METERS 


100 


26  Jun  1991 


SITE  fl  MflUI 


75 

50 

25 

0 

23 


TIME  (LOCRL) 


TIME  AVERAGED  CT^2  SITE  R  MRUI 


LOG  OF  CT^2 

Figure  34.  Echo-sounder  Turbulence  Profile, 26  June,  1991. 
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Figure  35,  Echo-sounder  Turbulence  Profile,  27  June,  1991 
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Figure  36.  Echo-sounder  Turbulence  Profile,  27  June,  1991. 
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Figure  39.  Echo-sounder  Turbulence  Profile,  27  June,  1991. 
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Figure  41.  Echo-sounder  Turbulence  Profile,  28  June,  1991. 
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Figure  42.  Echo-sounder  Turbulence  Profile,  28  June,  1991. 
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Figure  44.  Echo-sounder  Turbulence  Profile,  29  June,  1991. 
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Figure  55.  Echo-sounder  Turbulence  Profile,  01  July,  1991. 


97 


1  Jul  1991  SITE  fl  RMOS 


(f) 

OC 


u 

1- 

Li 

r 

z 


u 

o 

2 

(T 

OC 


TIME  (LOCAL) 


TIME  AVERAGED  CT-'-a  SITE  A  AMOS 


LOG  OF  CT--2 


Figure  56.  Echo-sounder  Turbulence  Profile,  01  July,  1991 
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Figure  57.  Echo-sounder  Turbulence  Profile,  01  July,  1991. 
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Figure  58.  Echo-sounder  Turbulence  Profile,  01  July,  1991. 
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Figure  59.  Echo-sounder  Turbulence  Profile,  02  July,  1991. 
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Figure  60.  Echo-sounder  Turbulence  Profile,  02  July,  1991. 
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Figure  61.  Echo-sounder  Turbulence  Profile,  02  July,  1991. 
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Figure  62.  Echo-sounder  Turbulence  Profile,  02  July,  1991 
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APPENDIX  C.  WIND  PROFILE  DATA  FOR  AMOS 


In  addition  to  Cy®  data,  wind  profile  dataware  also  collected  using  a  Doppler 
algorithm.  The  following  is  a  representative  sampling  of  this  data. 
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Figure  63.  Echo-sounder  Wind  Profile,  28  June,  1991. 
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gure  64.  Echo-sounder  Wind  Profile,  28  June,  1991. 


Figure  65.  Echo-sounder  Wind  Profile  29  June  1991. 
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gure  66.  Echo-sounder  Wind  Profile,  29  June,  1991. 
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Figure  67.  Echo-sounder  Wind  Profile,  29  June,  1991. 


107 


(vjl)  tySxB 


Figure  70.  Echo-sounder  Wind  Profile,  30  June,  1991. 


Figure  71.  Echo-sounder  Wind  Profile,  30  June,  1991. 
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gure  72.  Echo-sounder  Wind  Profile,  01  July,  1991. 


Figure  73.  Echo-sounder  Wind  Profile,  01  July,  1991. 
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igure  74.  Echo-sounder  Wind  Profile,  02  July,  1991. 


Figure  75.  Echo-sounder  Wind  Profile,  02  July,  1991. 


Ill 


LIST  OF  REFERENCES 

1.  Hall,  F.  F.  Jr.,  "Temperature  and  Wind  Structure  Studies  by  Acoustic  Echo- 
Sounding,"  Remote  Sensing  of  the  Troposphere,  U.  S.  Department  of 
Commerce  National  Oceanic  and  Atmospheric  Administration,  15  August, 
1970. 

2.  Neff,  W.  D.,  "Quantitative  Evaluation  of  Acoustic  Echoes  from  the  Planetary 
Boundary  Layer,"  A/OAA  Technical  Report  ERL  322-WPL  38,  June  1975. 

3.  Mousley,  T.  J.  and  Cole,  R.  S.,  "High  Frequency  Atmospheric  Acoustic 
Sounders,"  Atmospheric  Environment,  Vol.  13,  pp.  347-350,  1979. 

4.  Asimakopoulos,  D.  N.,  Helmis,  C.  G.  and  Stephanou,  G.  J.,  "Atmospheric 
Acoustic  Minisounder,"  Amer/can  Meteorological  Society,  pp.  345-350,  June 
1987. 

5.  Walters,  D.  L,  Private  Communication,  Naval  Postgraduate  School, 
Monterey,  California. 

6.  Moxcey,  L  R.,  Utilization  of  Dense  Packed  Planar  Acoustic  Echosounders 
to  Identify  Turbulence  Structures  in  the  Lowest  Levels  of  the  Atmosphere,  M. 
S.  Thesis,  Naval  Postgraduate  School,  Monterey,  California,  December  1987. 

7.  Weingartner,  F.  J.,  Development  of  an  Acoustic  Echosounder  for  Detection 
of  Lower  Level  Atmospheric  Turbulence,  M.  S.  Thesis,  Naval  Postgraduate 
School, Monterey,  California,  June  1987. 

8.  Wroblewski,  M.  R. ,  Development  of  a  Data  Analysis  System  for  the  Detection 
of  a  Lower  Level  Atmospheric  Turbulence  with  an  Acoustic  Echosounder,  M. 
S.  Thesis,  Naval  Postgraduate  School,  Monterey,  California,  June  1987. 

9.  McCrary,  J.  K. ,  High  Resolution  Cj^  and  Radial  Velocity  Measurements  Using 
a  High  Frequency  Monostatic  Acoustic  Echosounder,  M.  S.  Thesis,  Naval 
Postgraduate  School,  Monterey,  California,  June  1990. 

1 0.  Tatarski,  V.  I.,  Wave  Propagation  in  a  Turbulent  Medium,  Dover  Publications, 
New  York,  1961. 


112 


12.  Racine,  R.,  and  others,  "Mirror  Dome  and  Natural  Seeing  at  CFHT," 
Publications  of  the  Astronomical  Society  of  the  Pacific,  Vol.  103,  pp.  1020- 
1032,  September  1991, 

13.  Ochs,  G.  R.  and  Hill,  R.  J.,  "Optical  Scintillation  Method  of  Measuring 
Turbulence  Inner  Scale,"  Applied  Optics,  Vol.  24,  No.  15,  pp.  2430-2432,  1 
August  1985. 

14.  Hill,  R.  J.  and  Clifford,  S.  F.,  "Modified  Spectrum  of  Atmospheric 
Temperature  Fluctuations  and  its  Application  to  Optical  Propagation," 
Journal  of  the  Optical  Society  of  America,  Vol.  68,  No.  7,  pp.  892-899,  July 
1978. 

15.  Bass,  H.  E.,  "Atmospheric  Absorption  of  Sound;  Update,"  Journal  of  the 
Acoustic  Society  of  America,  Vol.  88,  No.4,  October  1990. 

16.  Olmstead,  M.  R.,  Development  of  a  Differential  Temperature  Probe  for  the 
Measurement  of  Atmospheric  Turbulence  at  All  Levels,  M.  S.  Thesis,  Naval 
Postgraduate  School,  Monterey,  California,  December  1988. 

17.  Fried,  D.  L,  "Optical  Resolution  Through  a  Randomly  Inhomogeneous 
Medium  for  Very  Long  and  Very  Short  Exposures,"  Journal  of  tthe  Optical 
Society  of  America,  Vol.  56,  No.  10,  pp.  1372-1379,  October  1966. 

18.  Oldenettal,  J.  R.,  "Note  on  the  Conversion  of  Seeing  Angle  to  r^  for  CIS 
Stellar  Image  Data,"  Private  Communication,  Air  Force  Maui  Observation 
Site,  Maui,  Hawaii,  4  September  1991. 

19.  Hoover,  C.  R.,  Investigation  of  a  Single  Point  Temperature  Probe  for 
Measurements  of  Atmospheric  Turbulence,  M.  S.  Thesis,  Naval  Postgraduate 
School,  Monterey,  California,  December  1991. 


113 


INITIAL  DISTRIBUTION  LIST 

No.  Copies 

1 .  Defense  Technical  Information  Center  2 

Cameron  Station 

Alexandria,  VA  22304-6145 

2.  Library,  Code  52  2 

Naval  Postgraduate  School 

Monterey,  CA  93943-5002 

3.  Richard  Frosch  1 

Code  PL/LIM 

Phillips  Lab 
Kirtland  AFB 
Albuquerque,  NM  87117 

4.  CPT  Ann  Slavin  1 

Code  PULTE 

Phillips  Lab 
Kirtland  AFB 
Albuquerque,  NM  87117 

5.  CPT  Andrew  Terzakis  1 

Code  PLyWE 

Phillips  Lab 
Kirtland  AFB 
Albuquerque,  NM  87117 

6.  Donald  L.  Walters  7 

Code  PH/WE 

Naval  Postgraduate  School 
Monterey,  CA  93943 

7.  LT  Timothy  S.  Mattingly  1 

246  Guthrie  Drive 

Bardstov'n,  KY  40004 


114 


